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ABSTRACT 

Aims. This study aims to determine the impact of stellar binary companions on the lifetime and evolution of circumstellar disks in 
the Chamaeleon I (Chal) star-forming region by measuring the frequency and strength of accretion and circumstellar dust signatures 
around the individual components of TTauri binary stars. 

Methods. We used high-angular resolution adaptive optics JHK^L' -bmid photometry and 1. 5-2.5 |^m spectroscopy of 19 visual binary 
and 7 triple stars in Chal - including one newly discovered tertiary component - with separations between ~25 and ~1000au. The 
data allowed us to infer stellar component masses and ages and, from the detection of near-infrared excess emission and the strength of 
Brackett-y emission, the presence of ongoing accretion and hot circumstellar dust of the individual stellar component of each binary. 
Results. Of all the stellar components in close binaries with separations of 25-100 au, 10^^^% show signs of accretion. This is 
less than half of the accretor fraction found in wider binaries, which itself appears significantly reduced (~44%) compared with 
previous measurements of single stars in Chal. Hot dust was found around 50^^^% of the target components, a value that is indis- 
tinguishable from that of Chal single stars. Only the closest binaries (< 25 au) were inferred to have a significantly reduced fraction 
(<25%) of components that harbor hot dust. Accretors were exclusively found in binary systems with unequal component masses 
Msecondary/Mprimary < 0.8, implying that the dctcctcd accelerated disk dispersal is a function of mass-ratio. This agrees with the finding 
that only one accreting secondary star was found, which is also the weakest accretor in the sample. 

Conclusions. The results imply that disk dispersal is more accelerated the stronger the dynamical disk truncation, i.e., the smaller the 
inferred radius of the disk. Nonetheless, the overall measured mass accretion rates appear to be independent of the cluster environment 
or the existence of stellar companions at any separation >25 au, because they agree well with observations from our previous binary 
study in the Orion Nebula cluster and with studies of single stars in these and other star-forming regions. 
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^T ■ 1, Introduction in terms of disk size, mass, and irradiation. Observational ex- 

^*^ ' ploration of these systems helps to advance the understanding 

rn Binary stars are among the most important branches of stellar of ^isk structure and evolution: the distribution of material in 

. . formation. Owing to their observationally complex nature, the the disk, dust particle dynamics, and grain growth leave an im- 

> ; impact of a stellar companion on the evolution of the circumstel- p^nt in the photometric and spectroscopic signatures of the star- 

lar disks, however, has yet to be explored in detail. Analytical ^isk systems. These can be compared with either predictions of 

and numerical studies consistently predict the outer truncation disk evolution mode ls (e.g., Kiev et al. 2008; Zsom et al. 2011; 

^ radius of a circumstellar disk to be less than one half of the iMtiller & Kle^l2Q12h or the distribution of planets in binary sy s- 

bin^y separation (.Artymowicz & Lubow||l994|; [Armitage et al.| tems (iRoell et al.ll2012) . Comparison with the latter constrains 






Il999|; iPichardo et al.| |2005|). Since typical binary separations pianet formation scenarios by reflecting whether the conditions 



^^n the order of a few ten to a few 100 au (|Leinert et al.| in the modified disks, e.g., evolutionary timescales, are favorable 

1 19931: [Mathieu 1994), disks around binary star components are fQj. ^^e formation of planets or not. 

significantly smaller than their 100-1000 au single star coun- Results from previous studies of visual TTauri binary stars 

terparts - a prediction qualitatively confirmed by the reduced ^^.^w an inconclusive picture: While the correlation of disk 

flux in submm binary observations ([Jensen et al. | ll994i [19961 ^i^e and mass with time until the disk disperses is expected 

[Harris etal.|!2012). The resultmg systems allow one to explore f^^^ evolutionary models (e.g., Monin eTaLl l2007h . observa- 

disk evolution with parameters that are distinct from single stars ti^ns appear to return no clear signature. For example, mixed 

binary systems consisting of one component harboring a disk 

* Based on observations collected at the European Organisation for and one bare stellar component exist in both configurations - 
Astronomical Research in the Southern Hemisphere, Chile. ESO Data the more massive or less massive component can be the disk- 
ID: 086.C-0762. bearing compone nt (White & Ghez 2001; Hartisan & KenvoiJ 
** e-mail: daemgen@astro . utoronto . ca 120031: IPrato et al.ll2003i:iMonin et al.ii2007i:iDaemgen et al.ii2012i: 
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ICorreia et al.ll2013l submitted). Accordingly, the larger and more 
massive disk or smaller and less massive disk may disappear 
first. As in the single-star case, the individual binary components 
appear to lose the ir inner disk between < 1 Myr and ~10Myr 
(iKraus et al.ll2Q12h . 

Nevertheless, binary disk lifetimes are typically shorter than 
in single stars and the less rnassive component loses its disk - on 
average - first (JMonin et al.ll2QQ7b . Most of the signatures, how- 
ever, are blurred or may even remain undiscovered since the typ- 
ically small sample sizes require combining data from diff'erent 
star- forming regions and from a wide mass range. This results in 
degeneracies when one tries to separate mutual influences in the 
parameter space, as stellar masses and the cluster environment 
are known to have significant influence on disk evolution. Larger 
sample sizes could so far only be investigated for spatially unre- 
solved systems that measure, e.g., the total disk mass in a binary 
system. While binary tren ds, suche as an accelerated disk dis - 
persal, could be confirmed (JKraus et al.ll2Q12trCieza et ani2009l) . 
variations with component mass and age as well as difl'erential 
component disk evolution were not accessible to these studies. 

Only few studies tackled these limitations by observing un- 
biased samples with coherent observational methods. Dedicated 
spatially resolved binary observations of the dense Orion Nebula 
cluster (ONC; Daemgen et al. 2012; Correia et al. 2013) and the 
young Taurus association (White & Ghez 2001) feature results 
which, however, require substantiation with data from other star- 
forming regions. On one hand, new observations can be used to 
decide which characteristics are universal and which are features 
of the individual regions under consideration. On the other hand, 
a careful selection of the target samples may allow one to extract 
information about disk properties as a function of crucial param- 
eters such as, e.g., age and mass. 

To expand the size of an available coherent dataset of disks 
in visual binaries, this paper focuses on the Chamaeleon I (Cha I) 
region. As part of the Chamaeleon complex, the Cha I dark cloud 
is close to the Sun (160±15pc; Whittet et al. 1997) with an age 
distribution peaking between 2 and 3 Myr (jLuhmanll2007b . Of 
the ~250 currently known stellar and subs tellar members ~50% 
show evidence of a circumstellar disk (JLuhman et al.l l2008l) . 
Close to 40 binaries and higher order multiples are known in a 
separation range between 16 and 1000 au, which implies a mul- 
tiplicity fra ction of about 30% - h igher by a factor of -1.9 than 
in the field (iLafreniere et al.ll2QQ8l) . Not only does Chal's prox- 
imity grant access to spatially resolved visual binaries with sepa- 
rations as close as a few tens of au through using adaptive optics 
on 8 m-class telescopes, but its older age compared with Taurus 
and the ONC can also provide first insights on disk evolution 
with time. 

The study presented in this paper explores a sample of 26 
multiple stars in the Chamaeleon I star- forming region. The pres- 
ence of accretion and/or hot dust in the circumstellar disks 
is evaluated for each individual stellar component using high- 
angular resolution near-infrared photometry and spectroscopy. 
This study complements our previous obs ervations of the ONC 
(iDaemgen et al.ll2012l: ICorreia et al.ll2013h to explore the evolu- 
tion of circumstellar disks in binaries as a function of age and 
environment. The observations are part of our ongoing eff'orts to 
characterize disk evolution as a function of stellar and binary pa- 
rameters, with the goal of composing a coherent picture of star 
and planet formation in multiple stellar systems. 



2. Observations and data reduction 

2.1. Sample definition 

Nineteen binaries and seven triple stars, including one 
newly discovered tertiary component to T33B, were ob- 
served in the ChamaeleonI star- forming region, with both 
SINFONI/VLT H-\-K-hsind integral field spectroscopy and 
close-in-time NACO/VLT JHK^L imaging between December 
2010 and April 2011. Twenty-f our of the targe t s were se- 
lected from the binary survey by 'L afreniere et al.l (12008) and 
two likely hosts of sub-stellar companions were additionally 
included (CHXR73, Luh man et al.1 l2006l: T 14. ISchmidt et al] 
■2008.) . This sample of 26 multiples is designed to contain nearly 
all multiples in Cha I in the separation range of 0'.'2 to ~6'' 
known when the targets were selected, which is equal to ~30- 
lOOOau at the distance of Cha I. The separation range was cho- 
sen to ensure high-quality spatially resolved observations of 
the individual binary components with ground-based adaptive 
optics-assisted spectroscopy. The brighter component of each 
binary was required to be brighter than y= 16.7, 7?= 17, and 
^= 13.5 to serve as the adaptive optics guide star for our NACO 
and SINFONI observations. All multiples were selected to have 
at least one component with a spectral type (SpT) in the T Tauri 
range, i.e., G-type or later. While typically all components of 
each multiple fulfill this criterion, two binaries have primary 
stars of comparably early spectral type (T 26 A, GO; T 41 A, B9). 
These components are excluded from some of the discussion in 
Sect. [3] to minimize the impact of biases caused by their rel- 
atively high mass. Likewise, the two substellar companions to 
CHXR 73 and T 14 were excluded because they turned out to be 
too faint for most of the spectroscopic evaluation. 

All multiples are confirmed members of Cha I, as inferred 
from their space velocity and indicators of youth (position in H- 
R diagram, infrared-excess, lithium abundance, stellar activity, 
iLuhm an 2004 ; iLuhman et al. 20 0l|). Physical binar ity was stud- 
ied by La freniere et al.l (2008) and JVogt et aP (12012.) through as- 
sessing the local stellar density and co-motion of binary com- 
ponents. While these authors found the probability of chance 
alignments to be < 10"^ for all but the widest binaries, physical 
association of the wide T 6, T 26, and T 39 multiples was con- 
firmed through multi-epoch imaging data. In addition, the spec- 
troscopy from this study allows us to rule out background gi- 
ants, and we did not detect significant changes in separation and 
position angle for any multiples of the sample when we com- 
pared with previously published measurements. The previously 
undetected companion to T33 B is discussed and shown to be 
physically related in Sect. 13.51 In the following, we will treat all 
multiples as physically bound systems. 

In every multiple, the brighter component in ^s-band is re- 
ferred to as the primary or 'A" component. This typically means 
that they are of earlier spectral type and thus higher mass than the 
other members of each system. In two cases (T 39, Hn4) the in- 
ferred mass of the B component is higher than that of the bright- 
est component in ^s- T39B is a spectrosc opic binary that we 
were unable to resolve (iNguyen et al.ll2012l) . Its mass is accord- 
ingly very likely overestimated and T 39 B not the most massive 
component of this higher-order multiple system. To be consis- 
tent with previous publications, we kept the identifiers as they 
are (the brighter component in K^ is referred to as A) but treated 
the more massive component as the primary for the derivation 
of, e.g., mass ratios^^ 

Nguyen et M] (l2012h examined all but five of our targets 
for spectroscopic binarity and found CHXR 28 B, T31 A, and 
T39B to show signs of spectroscopic binarity (see Table [TJ. 
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iNguyen et alJ also reported CHXR 47 and Hn 4 as spectroscopic 
binary candidates. These authors found, however, that the known 
visual companions are likely identical to the spectroscopic com- 
panions in their data. Since it is not possible to measure separate 
spectral or photometric properties of the spectroscopic compo- 
nents from our observations, the combined properties were used 
where applicable and possible implications of the binary nature 
are highlighted. These targets are excluded from those parts of 
the discussion where knowledge of the individual component 
properties is required. 

Five of the six triple stars in this sample show a clear hier- 
archical structure. Whenever binary parameters (such as separa- 
tion, mass ratios) are discussed, these triples are included as two 
binary stars, one consisting of the close binary system, the other 
one using the combined characteristics of the close binary (to- 
tal system mass, center position, etc.) to compose a binary with 
the wide tertiary component. One triple, T 39, does not show a 
distinct hierarchical structure and is either seen in a proj ection 
that veils an existing hierarchy, or its components are currently 
in an unstable non-hierarchical configuration (for a discussion 
see Lafreniere et al. 2008). We include this system as two bina- 
ries consisting of the AB and AC components, respectively (see 
Fig. [T] for the configuration). 

A list of the targets, spectroscopic binarity, and the dates of 
the observations can be found in Table [T] 

2.2. NACO imaging: pliotometry and astrometry 

2.2.1 . Observations and data reduction 

All targets were observed with NA OS-CONICA (NACO; 
iLenzen et alll2003l: iRousset et ani2003h imaging on UT4 of the 
Very Large Telescope (VLT). The adaptive optics system NAOS 
always used the brightest multiple component (or the binary it- 
self for close systems) as a guide star. All components of each 
target were observed in the same field of view in four broad-band 
filters (JHK^U) during a single night to minimize the impact 
of variability. Observations in JHK^ use the N90C10 dichroic 
and the S 13 camera with a pixel scale of 13.26mas/pix, U uses 
the IR wavefront sensor and the L27 camera with 27.1 mas/pix. 
Typical values of measured FWHM of our point sources are 
a.'074 in /, a'070 in H, a.'077 in ^s, and a.^l in L\ The dero- 
tator was used to keep the north/south direction aligned with 
the columns of the detector. All targets were observed in a five- 
position on-source dither pattern to allow for sky subtraction. 
Total integration times per target and filter are between 50 and 
800 sec in JH K^-hsind and equal to 176 sec in L\ with the ex- 
ception of the substellar companion host CHXR 73 of 528 sec. 
Photometric standard stars were always observed during the 
same night with the same instrumental setup as the science ob- 
servations as part of the standard NACO calibration plan. 

Pipeline-reduced data products were provided in the ESO 
data package for all science targets. A subsample was com- 
pared to images manually reduced with the NACO pipeline in 
Gasganc^ with the recommended settings to confirm their us- 
ability for photometric and astrometric analysis. The diff'erence 
between the two reductions was negligible (photometry, noise), 
thus the provided data products were used for the photometry. 
Manual pipeline reduction, however, was required for the pho- 
tometric standards, since no reduced data were available. For 
both the standards and science targets, the pipeline performs 
flat-fielding using lamp-flats and sky subtraction. After all the 



dithered images per target and filter were shifted and averaged, 
residual bad pixels were removed by substitution with the me- 
dian value of the surrounding good pixels. 

Images of all spatially resolved previously known multiples 
in the ^s-band are shown in Fig. [T] The close binary compan- 
ions in the triple systems T 26 and T 3 1 could not be resolved. 
Their binarity is revealed only by their elongated shapes (Fig.O 
in these observations with a FWHM of ~0'.'07. The estimated 
separations are <a.'05 (T26B) and ~a.'07 (T31B). The tertiary 
component to T 33 has not been reported before (Fig. E} and its 
physical association with the system is assessed in Sect. 13.51 

2.2.2. Photometry and astrometry 

Using a combination of calibrated and diff'erential component 
photometry, JHK^L photometry of all target components was 
measured. 

Diff'erential photometry was obtained through simultaneous 
aperture photometry of pairs of stars in the image (typically the 
two binary components). Triple stars allowed using accurate PSF 
photometry with the daophot task in IRAF using the PSF of 
the wide tertiary component as a reference. Despite possible iso- 
planatic eff'ects of the adaptive optics correction, subtracting of 
the flux- scaled reference resulted in residuals with typical am- 
plitudes ^1% of the original component flux. 

Separations and position angles, were inferred from the re- 
sults obtained with IRAF daof ind, which returns stellar posi- 
tions with an accuracy of a fraction of a pixel, or <0'.'01. The 
derotator was used to keep north up, so the inferred position an- 
gles are assumed to be absolute with a systematic uncertainty of 
0.5°. 

The photometric calibration was performed with photomet- 
ric standard stars observed during the same night or a night as 
close in time as possible, where the latter case is validated by 
the stablility of the zeropoint over the duration of several nights 
(rms error~0.02mag). Apertures of ~0'.'5-lV9 radius were used, 
either enclosing one binary component or, for close compan- 
ions, both components of a binary. We extracted the standard star 
photometry with the same individual apertures to avoid aperture 
corrections. Errors of the derived zeropoints caused by the vari- 
able adaptive optics PSF wings were estimated to be < 0.05 mag. 
Apparent magnitudes were then calculated considering normal- 
ized instrumental magnitudes, airmasses, Paranal extinction co- 
efficients, and photometric zeropoints. The resulting photometry 
is listed in Table [2] together with the measured separations and 
position angles. 

Two stars in the sample, T 26 A and T 33 A, turned out to be 
too bright for observations in L'-band with NACO. Their satu- 
rated L photometry cannot be used for further evaluation. The 
brown dwarf companion candidate T 14 B remains undetected in 
/. 

All JHK^ photometry presented in this section was ob- 
served in or converted to the 2MASS system ( Cutri et al. 2003), 
L p hotometry was converted to the ESO system when neces- 
sary (Ivan der Bliek et al.]|l996l) . Photometry with NACO needs 
no transformation to be compatible with this system and the 
photometry of the standard stars is available in these filter^. 
For all other necessary transformations of JHK^L photometry 
(e.g. CTTS locus, ext incti on vector; Sect. 13.11) the recipes of 
iBesseU & Bretj (Il988b and lCarpente^ (l200lh were used. 



v2.4.0, |http : //www . eso . org/sci/sof tware/gasganoT] 



|http : / /www . eso . org/s ci/facili ties/paranal/| 



|instrumerits/naco/tools/imaging_standards/nac~lw_lp . html | 





Table 1. Chal targets and observations 




















u 


RA 


DEC 








pre- 


telluric 


MJD 


phot. 


ref. 


MJD 


C/3 


rh m s\ 


C) 


Identifier 


Simbad Name^ 


SB^ 


optics^ 


reference 


spectroscopy 


JHK, 


L 


photometry 


_(D 


10 55 59.759 


-77 24 40.07 


T3 


AssChaT2-3 




0.25 


Hip 061 193 


55571.32071327 


S708-D_9141 


HD 75223 


55610.19226973 


a 


10 58 16.774 


-77 17 17.06 


T6 


Ass Cha T 2-6 


no 


0.25 


Hip 064501 


55559.34191219 


S372-S_9137 


HD 77281 


55602.28873422 


;s 


11 01 18.750 


-76 27 02.54 


CHXR9C 




no 


0.1 


Hip 061 193 


55571.34166417 


S372-S 


HD 77281 


55557.29178934 


c3 














Hip 050126 


55587.20283114 








H 
H 


11 02 32.654 


-77 29 12.98 


CHXR71 




no 


0.1 


Hip 050126 


55542.34920563 


HD 110621 


GL347A 


55641.08828209 














Hip 033708 


55565.35663229 










11 04 09.090 


-76 27 19.38 


T14 


Ass ChaT 2-14 


no 


0.25 


Hip 060429 


55590.21242836 


S708-D_9141 


HD 106965 


55590.34649026 


-§ 


11 05 14.673 


-77 11 29.09 


Hn4 




(SB2?)^ 


0.025 


Hip 040105 


55592.29434212 


HD 110621 


GL347A 


55641.10665852 


g 


11 05 43.002 


-77 26 51.75 


CHXR 15 






0.1 


Hip 050126 


55587.21137326 


HD 110621 


FS129 


55653.12026904 


1 


11 06 28.774 


-77 37 33.16 


CHXR73 






0.1 


Hip 050126 


55587.22603233 


HD 110621 


FS129 


55646.12072363 


^ 














Hip 067042 


55588.29378704 








11 07 20.722 


-77 38 07.29 


T26 


Ass Cha T 2-26 


no 


0.25 


Hip 060429 


55590.15206476 


HD 110621 


HD 106965 


55653.14097253 




11 07 28.256 


-76 52 11.90 


T27 


Ass Cha T 2-27 




0.1 


Hip 067042 


55678.05304754 


HD 110621 


HD 106965 


55653.16332159 


1 


11 07 55.887 


-77 27 25.75 


CHXR 28 




CHXR28B=SB1? 


0.1 


Hip 064501 


55594.23506763 


HD 110621 


HD 106965 


55653.18627437 


11 08 01.486 


-77 42 28.85 


T31 


Ass ChaT 2-31 


T31A=SB2? 


0.1 


Hip 065181 


55569.28001642 


HD 110621 


HD 106965 


55653.20674135 


•Ji 


11 08 02.975 


-77 38 42.59 


ISO 126 


ISO-Chal 126 


no 


0.1 


Hip 040105 


55589.31037568 


HD 110621 


HD 106965 


55653.22875679 




11 08 14.938 


-77 33 52.19 


T33 


Ass Cha T 2-33 


no 


0.25 


Hip 060429 


55590.16543515 


HD 60778 


HD 60778 


55666.03397411 


i 


11 09 11.723 


-77 29 12.48 


T39 


Ass Cha T 2-39 


T39B=SB1 


0.25 


Hip 060429 


55590.20255577 


HD 60778 


HD 60778 


55666.05782756 


^ 


11 09 18.129 


-76 30 29.25 


CHXR 79 




no 


0.1 


Hip 061 193 


55589.34228160 


HD 60778 


HD 60778 


55666.08169563 


^ 


11 09 50.019 


-76 36 47.72 


T41 


Ass ChaT 2-41 




0.1 


Hip 067969 


55589.26156993 


HD 110621 


HD 75223 


55605.33691345 


^ 


11 09 54.076 


-76 29 25.31 


T43 


Ass Cha T 2-43 


no 


0.1 


Hip 040105 


55589.33528044 


HD 60778 


HD 60778 


55666.11141484 


s 


11 09 58.738 


-77 37 08.88 


T45 


Ass Cha T 2-45 


no 


0.1 


Hip 067969 


55589.27490141 


HD 60778 


HD 60778 


55666.13248263 


t 


11 10 38.018 


-77 32 39.90 


CHXR 47 




(SB2)^ 


0.025 


Hip 040105 


55592.28742058 


HD 60778 


HD 60778 


55666.15952256 


11 11 54.002 


-76 19 31.14 


CHXR49NE 






0.1 


Hip 061 193 


55589.36918014 


HD 110621 


HD 75223 


55606.24352530 


^ 


11 12 24.415 


-76 37 06.41 


T51 


Ass ChaT 2-51 


no 


0.1 


Hip 067969 


55589.29062811 


HD 60778 


HD 60778 


55670.07307583 




11 12 42.689 


-77 22 23.05 


T54 


Ass Cha T 2-54 


no 


0.1 


Hip 040105 


55589.32408822 


S372-S_9137 


HD 77281 


55602.31526911 


'^ 


11 14 25.890 


-77 33 06.10 


Hn21 




no 


0.25 


Hip 040105 


55592.32103162 


HD 110621 


HD 75223 


55606.29302649 




11 14 50.319 


-77 33 39.04 


B53 


BYB53 


no 


0.1 


Hip 040105 


55589.30224627 


HD 110621 


HD 75223 


55606.34424266 




11 18 20.600 


-76 21 58.00 


CHXR 68 




no 


0.25 


Hip 060429 


55590.19058466 


HD 60778 


HD 60778 


55666.18125080 


e 
^ 


Notes. ^^^ Identifier to be used with [ 
^'^ SINFONI pre-optics used (see text 








diff'erent from identifier in the 3rd column. ^^^ Spectroscopic binarity referenc 
and CHXR 47 are probably identical to the visual companions in this paper (s 




a 


Mmbad ( http : / 


/simbad . u- strasbg . fr/simbad/ ) if 


e: 'Nguven et al, 


(I2OI2I1 


iri 


for definition) *^^ 


The spectroscopic companions to Hn 4 


ee.Nsuyen et al. 


20T2I). 
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B53 

« 

• 

1" 


CHXR 9C 

Bb 

Ba 'a 


CHXR 15 

• 

• 


CHXR 28 

Ab 

i 

Aa 

• 
B 










CHXR 47 


CHXR49NE 
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A 


CHXR 68 

Ab. 
Aa 

B 


CHXR 71 

• 








chxrv'b 

B - , 


CHXR 79 

* 


Hn 4 

B 
• • 
A 


Hn 21 

• 








ISO 126 


T3 

m 

m 


T6 


T14 








T27 

• 


T39 

c 

B 

m 

A 


T41 

• 


T43 

• 








T45 


T51 

m 


T54 













Fig. 1. NACO ^s-band imaging of all target binaries and triples (except T26, T31, and T33 in Fig.O. Components are marked 
with their component identifier when it is not clear whether they are primary or secondary. The bar in each subpanel has a length of 
V\ The color stretch is linear. North is up and east to the left. 
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Fig. 2. Close binary components of T 26 (left) and T 31 (middle), i.e., the outhem component in both panels, could not be spatially 
resolved with NACO ^s-band imaging despite the good FWHM of ~(y.'070. From the apparent elongation we estimate that the 
separation of T26BaBb is not much larger than (/.'OS and just at the FWHM for T 31 BaBb. The ^g-band image of T33 (right) 
shows t he n ewly discovered faint tertiary component close to T 33 B. The physical association of all three components is discussed 
in Sect. 13.51 Image stretches and orientation are the same as in Fig.[TJ 



Some targets were observed multiple times for various rea- 
sons like bad weather, poor Strehl correction, or a missing pho- 
tometric reference in one or more filters. The subset of repeti- 
tions that returned data of good quality in one or more filters 
were reduced and measured as described. More than half of all 
targets with repeated measurements exhibited significant pho- 
tometric variability of up to -0.2 mag and diff'erential magni- 
tude diff'erences of <0.1mag. A scatter of this magnitude can 
be caused by several mechanisms, including intrinsic variability 
through chromospheric activity, variable accretion, or occulta- 
tion from a warped in ner disk as a result of magneto spheric star- 
disk interaction (e.g ., ICarpenter et al.ll2QQlUBouvier et aO 19991: 



disk mteracti on (e.g . , 
lAlencaretall2QlQh . 



Despite occasional re-observations in individual filters, for 
all targets only one complete set of observations exists, i.e., all 
of JHK^L were observed during one night together with photo- 
metric references. These observations were used for subsequent 
evaluation. 



2.3. SINFONI H+K-band IFU spectroscopy 

2.3.1 . Observations and data reduction 

Integral field spectroscopy was obtain ed with SINFONI 
(Eisenhauer et al.ll2QQ3l:lBonnet et al.ll2QQ4 on the VLT. All ob- 
servations used the adaptive optics module, locked on the bright- 
est multiple component or close binary. The off'ered H-\-K con- 
figuration was used to obtain simultaneous spectroscopy from 
-1.5 to ~2.5 |Ltm for an array of 32x64 spatial pixels. Depending 
on the separation of the target binary, one of the 0.25, 0.1, or 
0.025 pre-optics was used, resulting in a final spatial pixel scale 
in the reduced images (after resampling) of OV 125/pix, O'.'05/pix, 
and O'.'0125/pix. The spectral resolution is R ^ 1500 with 
~5 A/pix in dispersion direction. The typical FWHM is -O'.'IO- 
0'.'25, sufficient to spatially resolve all binaries and wide tertiary 
components in the sample. All close components in triples, how- 
ever, are only marginally resolved and only the combined spectra 
were used for further evaluation. Except for CHXR73, all tar- 
gets were observed in a four-position dither pattern. In the two 
larger scales (0.25 and 0.1 pre-optics) all four dithers were on 
source; the two closest binaries (Hn4 and CHXR47) were ob- 



served with alternating on- and off'-source dithers with the 0.025 
pre-optics. The total integration time was 2 min per target. To 
increase the S/N of the faint secondary around CHXR73, this 
target was observed with 16 dither positions totaling 32 min of 
integration time. Each science observation was accompanied by 
the observation of a bright telluric standard star of known spec- 
tral type (B2 to B6) with identical system settings at similar air- 
mass. 

As suggested in the SINFONI data reduction cookboolu, we 
removed a pattern of dark lines from all science and calibration 
files - caused by a hard-coded bias-removal of the SINFONI in- 
strument - with the provided IDL script. The subsequent reduc- 
tion of all science and telluric reference frames used the standard 
SINFONI pipeline in EsoRe^ The pipeline executes bad-pixel 
removal, flat-fielding using a lamp flat, background subtraction, 
and wavelength calibration. A correction for difl'erential atmo- 
spheric refraction was applied to exposures using the pre-optics 
0.1 and 0.025. For each set of 4 or 16 dither exposures, the 
pipeline outputs a data cube that consists of 2172 fully reduced 
and combined 2D images, each corresponding to a wavelength 
between 14 750 and 24 930A, in steps of 5A. 

Before extracting the stellar ID spectra from the data cubes, 
unusable pixel data (residual bad pixels, nearby "negative stars" 
from the sky subtraction, or blending with other stars in the im- 
age) were replaced by the median value computed in a two pixel- 
wide annulus around the defect at the wavelength they occur. 
The accuracy of this method was tested by applying the median 
substitution to good regions in a test cube. The difl'erence be- 
tween the extracted flux from the regular and the partly substi- 
tuted cube was lower than 1% at all wavelengths. 

All spectra were extracted twice with slightly difl'erent algo- 
rithms for difl'erent purposes: (a) to be compared with template 
spectra, which required a good recovery of the overall spectral 
shape, and (b) to measure Bry emission with a good signal-to- 
noise ratio (S/N). 



^ |h ttp://www.eso.org/s ci/facilities/paranal/| 
|inst ruments/sinfoni/doc/1 
"^ |http : //www . eso . org/scT/sof tware/cpl/esorex . html] 
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Table 2. Individual component apparent magnitudes 





separation^ 


PA^ 




J' 


H^ 


K.' 


L'^ 


4 phot 


Ek.-l' 


Name 


[arc sec] 


[au] 


n 




[mag] 


[mag] 


[mag] 


[mag] 


[mag] 


[mag] 


B53 


0.28 


45 


235.8 


A 


10.80 + 0.03 


10.03 + 0.04 


9.87 ±0.02 


9.66 ±0.01 


0.4 ±0.5 


0.03 










B 


11.96 + 0.05 


11.44 + 0.06 


11.30 ±0.04 


10.83 ±0.03 


0.0 ±0.5 


0.21 


CHXR9C 


0.90 


144 


79.7 


A 


10.83 + 0.04 


10.01+0.03 


9.84 ±0.06 


9.76 ±0.02 


0.7 ±0.6 


-0.11 




0.12 


19 


69.5 


Ba 


11.15 + 0.04 


10.33 + 0.03 


10.16 ±0.06 


10.08 ±0.02 


0.7 ±0.6 


-0.12 










Bb 


11.82 + 0.04 


11.06 + 0.03 


10.84 ±0.06 


10.65 ±0.02 


0.5 ±0.6 


0.00 


CHXR15 


0.31 


50 


198.2 


A 


11.91+0.03 


11.30 + 0.04 


10.99 ±0.02 


10.40 ±0.04 


1.0 ±0.3 


0.20 










B 


12.38 + 0.04 


11.80 + 0.04 


11.53 ±0.02 


11.02 ±0.05 


1.1 ±0.4 


0.12 


CHXR28 


1.84 


294 


117.1 


Aa 


9.99 + 0.03 


9.03 + 0.04 


8.77 ±0.02 


8.62 ±0.02 


2.0 ±0.5 


-0.05 




0.12 


20 


7.6 


Ab 


10.44 + 0.03 


9.38 + 0.04 


9.09 ±0.02 


8.81 ±0.02 


2.7 ±0.5 


0.04 










B 


10.46 + 0.03 


9.44 + 0.04 


9.15 ±0.02 


9.00 ±0.02 


2.5 ±0.5 


0.17 


CHXR47 


0.17 


27 


328.9 


A 


10.34 + 0.03 


9.42 + 0.03 


8.87 ±0.03 


8.27 ±0.03 


3.6 ±0.2 


0.32 










B 


10.80 + 0.04 


9.83 + 0.03 


9.35 ±0.04 


8.77 ±0.03 


3.6 ±0.3 


0.29 


CHXR49NE 


0.27 


43 


48.6 


A 


10.99 + 0.03 


10.22 + 0.04 


10.02 ±0.02 


9.64 ±0.01 


1.4 ±0.3 


0.12 










B 


11.09 + 0.03 


10.37 + 0.04 


10.15 ±0.02 


9.73 ±0.01 


0.9 ±0.3 


0.14 


CHXR68 


4.40 


704 


212.1 


Aa 


10.59 + 0.03 


9.89 + 0.04 


9.61 ±0.02 


9.67 ±0.04 


0.4 ±0.5 


-0.22 




0.09 


15 


21.9 


Ab 


10.74 + 0.03 


10.07 + 0.04 


9.84 ±0.02 


9.71 ±0.04 


0.0 ±0.5 


-0.01 










B 


11.41 + 0.03 


10.67 + 0.03 


10.40 ±0.03 


10.28 ±0.03 


0.6 ±0.4 


-0.06 


CHXR71 


0.56 


90 


68.5 


A 


11.55 + 0.03 


10.69 + 0.04 


10.45 ±0.02 


10.24 ±0.02 


1.2 ±0.5 


-0.02 










B 


13.12 + 0.03 


12.41+0.04 


12.06 ±0.02 


11.27 ±0.03 


0.0 ±0.3 


0.53 


CHXR73 


1.29 


206 


234.6 


A 
B 


12.64 + 0.35 
17.84 + 0.37 


11.40 + 0.04 
16.35 + 0.10 


10.94 ±0.02 
15.47 ±0.16 


{...y 


4.5 ±2.5 
6.7 ±3.1 




CHXR79 


0.88 


140 


211.0 


A 


11.84 + 0.04 


10.31+0.03 


9.32 ±0.03 


8.03 ±0.03 


6.9 ±0.3 


0.76 










B 


13.84 + 0.05 


12.50 + 0.04 


11.80 ±0.04 


10.87 ±0.08 


5.6 ±0.4 


0.38 


Hn4 


0.17 


27 


295.2 


A 


11.65 + 0.03 


10.76 + 0.04 


10.45 ±0.02 


10.17±0.01 


1.7 ±0.5 


-0.07 










B 


11.70 + 0.03 


10.77 ±0.04 


10.47 ± 0.02 


10.32 ±0.01 


1.9 ±0.5 


-0.17 


Hn21 


5.48 


876 


68.5 


A^ 


12.09 + 0.03 


11.14 + 0.04 


10.70 ±0.02 


10.18 ±0.02 


3.5 ±0.3 


0.10 










B^ 


12.90 + 0.04 


12.09 ±0.04 


11.56 ±0.03 


10.77 ±0.03 


1.5 ±0.3 


0.38 


ISO 126 


0.28 


45 


228.8 


A 


12.55 + 0.03 


10.48 ±0.04 


8.87 ±0.02 


6.83 ±0.02 


10.5 ±0.3 


1.87 










B 


12.70 + 0.03 


11.21 ±0.05 


10.11 ±0.04 


8.33 ±0.02 


4.3 ±0.2 


1.41 


T3 


2.21 


354 


288.4 


A 


10.97 + 0.03 


9.64 ±0.02 


8.58 ±0.03 


7.12±0.01 


4.9 ±0.2 


1.05 










B 


11.62 + 0.03 


10.86 ±0.02 


10.33 ±0.03 


9.55 ±0.02 


1.2 ±0.5 


0.50 


T6 


5.12 


819 


122.3 


A 


8.76 + 0.04 


8.23 ±0.03 


7.70 ±0.05 


6.90 ±0.01 


0.0 ±0.5 


0.76 










B 


12.00 + 0.08 


11.65 ±0.06 


10.95 ±0.08 


10.39 ±0.02 


0.8 ±0.3 


0.26 


T14 


2.67 


427 


299.6 


A 
B 


9.45 + 0.02 


8.88 ±0.02 
15.01 ±0.36 


8.48 ±0.03 
16.54±4.16 


7.66 ±0.01 
13.01 ±0.15 


0.0 ±0.5 


0.72 


T26 


4.57 


731 


201.3 


A 


7.89 + 0.03 


7.06 ±0.04 


6.45 ±0.02 


sat 


0.0 ±0.5 












BaBb^ 


11.55 + 0.04 


10.75 ±0.04 


10.24 ±0.03 


9.90 ±0.04 


1.8 ±0.5 


0.01 


T27 


0.78 


125 


13.3 


A 


11.27 + 0.03 


10.48 ±0.04 


10.08 ±0.02 


9.38 ±0.03 


0.9 ±0.3 


0.50 










B 


11.68 + 0.03 


10.99 ±0.04 


10.77 ±0.02 


10.42 ±0.04 


1.1±0.3 


0.08 


T31 


0.66 


106 


182.6 


A 


9.32 + 0.03 


8.17 ±0.04 


7.38 ±0.02 


6.14 ±0.02 


3.0 ±0.3 


0.98 










BaBb^ 


9.74 + 0.03 


8.83 ±0.04 


8.51 ±0.02 


8.24 ±0.04 


1.8 ±0.5 


0.04 


T33 


2.40 


384 


284.4 


A 


9.65 + 0.03 


8.01 ±0.03 


6.68 ±0.02 


sat 


4.3 ±0.3 






0.11 


18 


85.6 


Ba 


10.09 + 0.03 


9.27 ±0.03 


9.07 ±0.04 


(.../ 


0.8 ±0.4 












Bb 


11.38 + 0.04 


10.56 ±0.03 


10.32 ±0.07 


(.../ 


1.0 ±0.6 




T39 


1.25^ 


200 


18.4 


A 


10.37 + 0.03 


9.67 ±0.03 


9.48 ±0.03 


9.88 ±0.03 


0.0 ±0.4 


-0.59 




4.50^ 


720 


69.9 


B 


10.77 + 0.03 


10.00 ±0.03 


9.83 ±0.03 


10.31 ±0.03 


0.4 ±0.4 


-0.66 










C 


11.01+0.04 


10.29 ±0.03 


10.09 ±0.04 


9.75 ±0.03 


0.2 ±0.5 


0.14 


T41 


0.78 


125 


324.3 


A 


7.61+0.03 


7.35 ±0.04 


7.28 ±0.02 


7.36 ±0.01 


0.0 ±0.3 


-0.03 










B 


11.25 + 0.28 


10.57 ±0.13 


10.27 ±0.06 


9.93 ±0.04 


0.3 ±2.5 


0.09 


T43 


0.78 


126 


352.0 


A 


11.56 ±0.03 


10.32 ±0.03 


9.55 ±0.02 


8.64 ±0.03 


5.5 ±0.3 


0.45 










B 


12.85 + 0.03 


11.75 ±0.04 


11.07 ±0.02 


10.11 ±0.03 


3.6 ±0.3 


0.51 


T45 


0.74 


119 


51.8 


A 


10.01+0.03 


8.93 ±0.03 


8.16 ±0.03 


7.06 ±0.03 


3.2 ±0.2 


0.78 










B 


11.78 + 0.03 


10.99 ±0.03 


10.59 ±0.03 


10.10 ±0.04 


2.3 ±0.2 


0.11 


T51 


1.97 


316 


161.8 


A 


9.32 + 0.03 


8.67 ±0.03 


8.28 ±0.02 


7.15 ±0.03 


0.0 ±0.5 


1.03 










B 


11.39 + 0.03 


10.72 ±0.03 


10.38 ±0.03 


9.66 ±0.03 


0.0 ±0.5 


0.53 


T54 


0.24 


38 


248.0 


A 


8.65 + 0.04 


8.24 ±0.03 


8.24 ±0.06 


8.19±0.01 


0.0 ±0.5 


0.02 










B 


10.16 + 0.05 


9.66 ±0.03 


9.65 ±0.06 


9.51 ±0.02 


0.0 ±0.5 


0.04 




uncertainties are 


Asep^a.'Ol, APA 


= 0.5°. Projected separations 


in au assume a 


. distance to Chal of 160pc 




Notes. ^"-^ Separation 


dWhittetetaLl 



1997). For hierarchical triples, separations between a wide component and the close pair are given as well as the separation of the close components. 
^^^ Uncertainties are statistical errors from the differential photometry. Systematic uncertainties from relative aperture photometry are ~0.1 mag. 
^^^ Due to a saturation of the reference star, no calibrated L photometry could be measured. The magnitude difference is AL' = L^-L^ = 4.69 + 0.19 
{d) jj^gsg components are also known as Hn21 E (primary) and Hn21 W (secondary). ^^^ These previously known visual binary components could 
not be resolved with our NACO adaptive optics observations. System photometry is listed. ^^^ These components are not separated in the L'-band 
images. The combined photometry returns Lg^^g^ = 8.60 + 0.03. ^^^ As no clear hierarchy is visible, sep. and PA of the secondary and tertiary 
components with respect to the primary are given. 
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(a) Extraction with preservation of the spectral slope. 
Because adaptive optics correction produces a wavelength- 
dependent PSF, we applied a custom extraction routine. To ac- 
commodate the curved trace of each star throughout the cube 
(the position of the star between 1.5 and 2.4 yum can differ by 
up to more than 1 spatial pixels) and the changing shape of 
the PSF as a function of wavelength, the center of each stel- 
lar component (target and standard) was traced and fitted by a 
fourth-order polynomial and its wavelength-dependent FWHM^ 
measured at every A by fitting a two-dimensional Moffat profile 
to the stellar disk. Multiple extractions were then conducted for 
each target component with aperture radii equal to 1, 2, 4, 6, 8, 
and 10 times the measured FWHM^. If the target PSF appeared 
to be strongly elliptical, or if two very close components of a 
multiple were extracted together, elliptical apertures were used, 
designed from an elliptical 2D Moffat profile fit. To make sure 
that the overall shape of the extracted spectrum was conserved, 
all extractions were compared the result from the largest aper- 
ture (lOxFWHM^), which contains close to 100% of the flux. 
All extractions that agreed with this reference spectrum to better 
than 5% at all wavelengths were then examined for their noise. 
The extraction with the highest S/N was selected (typically 2- 
4xFWHM^). 

(b) Extraction with maximum signal to noise. To determine 
the Bry flux, only local correction is necessary. Thus, the overall 
spectral shape does not need to be conserved and the extraction 
can be optimized for S/N in the Bry region. This was achieved by 
using small extraction apertures (0.5-1 xFWHM^). The typical 
improvement in S/N over the continuum-conserving extraction 
at Abvy is a factor of ~ 2. In the following, we refer to the results 
from (a) as flux-conserved extraction and from (b) as high-S/N 
extraction. 

Standard- star spectra were extracted with the same proce- 
dures. To produce pure telluric spectra, convolved only with the 
instrumental profile, intrinsic absorption lines (except Bry) were 
removed by interpolating the continuum. Since a careful recov- 
ery of Bry flux is crucial for measuring accretion, this absorp- 
tion line was re moved according to the procedure described in 
iDaemgen et al.l ^2012), which uses an estimate of the telluric 
spectrum around Bry from target spectra with neither absorp- 
tion nor emission in this line. We corrected for the intrinsic 
continuum shape of the standard stars by dividing blackbody 
spectra with Teff according to their spectral type. Telluric wave- 
lenghts were fine-tuned by applying subpixel shifts after cross- 
correlation of 20-pixel-wide sections of the telluric and target 
spectra. 

For two of our target multiples, CHXR9C and CHXR71, 
two epochs of data with good S/N were taken (see Table[T]). Both 
data sets were reduced and extracted according to the scheme de- 
tailed above. Since no significant variability was apparent in the 
data, both extractions of all respective components were aver- 
aged to increase the S/N. 

The final S/N of the reduced spectra depends on source 
brightness but never drops below 20 anywhere in the spectral 
windows of H and ^-band for the shape-conserved reductions 
except for the two substellar companions. The typical S/N per 
pixel around Bry in the high-S/N reduction is around 80. The 
spectral region of low atmospheric transparency between 18000 
and 20000 A has a low S/N and was not used in the following. 
The final reduced and extracted spectra of all multiple star com- 
ponents are shown in Appendix [A] the most important stellar 
absorption features (see Table [3]) and Bry are marked. 



Table 3. Spectral features identified in the observed spectra 



^ 


Width 






[A] 


[A] 


Species 


Transition 


14 877.5 




Mgl 


3^^^1,2,3-4^^23,4 


15 025.0 




Mgl 


As^S'l-Ap^Pl 


15 040.2 


* 


Mgl 


As^S\-Ap^Pl 


15 168.4 




KI 


M^D,i2-AfFli, 


15 740.7 




Mgl 


Ap^P^-Ad^Di 


15 765.8 


^ 


Mgl 


Ap^P^-Ad^D^ 


15 888.4 




Si I 


As^P^-Ap^Pi 


16719.0 




All 


'Pu2-'D3/2 


16750.5 


* 


All 


'^?/2-'^5/2 


17110.1 




Mgl 


As^So-Ap^P'i 



Notes. Characteristic photospheric features of late-type stars be- 
tween 14500 A and 18000 A. Features between 20000 and 25000 A 
are listed in Table 2 of [Daemgen et al. (2012). Asterisks mark 
lines that blend with the line with the next shortest wave- 
length (the previous entry in the list). Transitions selected from 
Gushing etal. (2005) and the NIST Atomic Spectra Database 
( http : //vivivi . ni St . go v/pml/data/asd . cfmj . 



2.3.2. Flux calibration 

After extraction, all spectra wer e flux-calibrated with a routine 
similar to the one described in Daemgen et al. (2012). ^s-band 
photometry was compared with synthetic photometry derived 
from convolving the uncalibrated SINFONI H-\-K spectra with 
a ^s filter curve, which is superior to an extraction with H or 
H-\-K filters in terms of noise and Bry calibration. Since the to- 
tal flux of the high-S/N spectra is not conserved, these use the 
calibration factor derived from the flux-conserved extraction of 
the same target. 



3. Results 

3.1. Color-color diagram: photometric extinctions and color 
excesses 

Fig. [3] shows the (J-H)-(H-K,) and (H-K,)-(K,-U) color- 
color diagrams for Chal binaries with photometric data from 
Table [21 The data are compared with the sequence of dwarfs and 
giants, the CTTS locus, and to a sequence of K4-M9 pre-main 
sequence stars taken from Luhman et al.l (I2Q1Q ). The data show a 
clear clustering around the tip of the PMS sequence at a spectral 
type -Ml. 

The color-color diagrams were used to derive photometric 
extinctions. The extinction Ay °^ of a target component was de- 
termined by moving it toward bluer colors along the reddening 
vector in the (J-H)-(H-K^) diagram until intersecting either 
the CTTS locus or the PMS sequence. To decide whether to 
deredden to the CTTS locus or the PMS sequence, the (H-Ks)- 
(K^-L') diagram was used since here the CTTS locus is clearly 
separate from the PMS, dwarf, and giant sequences. Targets that 
are clearly associated with the CTTS locus in {H-K^)-{K^-L), 
i.e., fall along it or cannot otherwise be dereddened to the PMS 
sequence, were dereddened to the corresponding CTTS locus in 
(J-H)-(H-Ks). All other target components were dereddened 
to the PMS sequence. Any targets with J-H <0.5 were as- 
signed an extinction of Av=0. Uncertainties were determined 
from dereddening according to the uncertainties of the color 
measurements. The photometric extinctions Ay °^ are listed in 
Tabled 
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Fig. 3. Color-color dia grams with photometric data from Table |2] Data are shown as red crosses, the giant and dwarf sequences (thin 
continuous lines) fr om lBessell & Bretj (Il988l) . the CTTS locus (iMeyer et al.ll 19971) . and a reddening vector of length Ay = 10 mag 
("Cohenet al."l98l') are overplotted. The dotted locus just above the dwarf sequence is the sequence derived from PMS stars in 
iLuhm an et al. (2 010 ). Targets between the two dash-dotted lines can be dereddened to the PMS locus. 



Component extinction s of s ome of the wider binaries 
are published in iLuhmanI (l2QQ7h . and values typically agree 
to within AAy^ +lmag (after conversion Aj = 0.276Ay; 
ICohen et al.lfl98lh . 

Color excesses were derived for all targets according to 



Ek,-l' = (^s - L'Us - 0.04Av - (K, - L% 



(1) ^ 



with referenc e (Kf^-L')o dwar f colors , converted to the 2MASS 
system, from Bessell & Bretj (1988) according to the spectral 
types in Table |4] (see Sect. 13.21) . It is possible that the derived 
excesses over- or underestimate the real excess, since the near- 
IR colors of low-mass PMS stars and dwarfs of the same spec- 
tral type differ by u p to a few times .1 mag (compare e.g. 
iBessell & Bret3 1 19881 with iLuhman et al.ll201 Q). Since there are 
no well-calibrated K^-L values for pre-main sequence stars 
available in the literature, we estim ated the size of this effect by 
comparing PMS (^-[3.6]) colors (Luhman et al.ll2QlQ|) to dwarf 
K^-L colors and found a difference between 0.05 and 0.17 mag 
in the spectral type range between K4 and M6. Although the 
Spitzer/IRAC [3. 6] -band and L are not identical, this gives an in- 
dication of the magnitude of this effect. We assigned systematic 
uncertainties of 0.2 mag to the derived Ek.-u values in Tabled 

3.2. Spectral types, extinction and veiling 

To measure spectral types, veiling, and to obtain another inde- 
pendent estimate of extinction, we compared the extracted flux- 
conserved spectra with dwarf spectra from the IRTF library in 
a simi lar way as described for ONC stars in iDaemgen et al.l 
(l2Q12h . Since Chal is older than the ONC, the spectra in this 
sample will be even more similar to dwarf spectra in the IRTF 
library. Indeed, the strong concentration of Cha I targets along 
the dwarf sequence in a iy(NaI-hCaI)-iy(CO[2-0]-hCO[4-2]) 
diagram (Fig. |4]) indicates the similarity between the depth of 
the photospheric absorption features and thus the surface grav- 
ity logg between dwarfs and the Chal PMS stars. Accordingly, 
dwarf spectra were assumed to be suited to serve as templates 
and were used to derive spectral parameters. 
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Fig. 4. Comparison of the Cha I target components with dwarfs 
(red continuous) and giants (blue dashed) from the IRTF spectral 
library in terms of their equivalent widths in the mostly tem- 
perature sensitive Nal and Cal lines and the surface-gravity- 
sensitive ^^CO band heads. The dwarf and gi ant loci are de- 
rived from the IRTF spectral library (see Fig. 4 in lDaemgen et al.l 
2012). Strong veiling reduces the measured equivalent widths 
and explains the clustering of targets close to the origin and to 
the left of the dwarf sequence. 



To limit the number of templates that need to be tested 
for their compatibility with each target spectrum, we ini- 
tially selected the spect ral types following the extinction- 
i ndependent Wilkin^ 's Q (IWilking et al.ll 19991) and IH2O indices 
(Comeron et al. '2000"). Since our observations do not cover 
the whole wavelength range used by the original definition of 
Q = Fl/F2(F3/F2y-^^ [with Fl (2.07-2.13 ^tm), F2 (2.276- 
2.285 |Lim), and F3 (2.4-2.5 |Ltm)], we redefined the last bin to 
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Fig . 5. Wilking's g and /h .o indices after I Wilkin g et al. I (Il999h 
and lComeron et al.' ('2000), measured for the dwarf templ ate se- 
quenc e from the IRTF library (Gushing et al. 2005; Rayner et al.l 
l2009l) . The dashed and dotted lines are best second-order poly- 
nomial fits to the data. 
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Fig. 6. Same as Fig.|4l but target components with a continuum 
significantly diff'erent from that of a dwarf star between 1 .45 and 
1.7 jLim are marked with red diamonds. Low- veiling components 
with vk < 0.1 are marked with black squares and high- veiling 
components r^ > 0. 1 with error bars only. 



be F3* =2.4-2.43 |LLm. The introduced wavelength dependence 
of Q is not significant in the limited spectral type range of the 
target sample, as was tested on template spectra of the same 
spectral types. A second-order polynomial was fit to both indices 
(Fig. [5]) derived from the templates as a fun ction of spectral type 
assigned by the authors of the sequence (ICushing et al.l 120051: 
iRayner et al.M2009.) . The two derived relations for subclasses of 
Mare 



with refined grids in extinction Ay, continuum excess k, and con- 
stant c according to 



(SpTM)e = 44.3-85.7 2 + 41.9 2^ , 
(SpTm)/,^, = -26.3 + 33.9/^,0-6.18 4^ 



(2) 
(3) 



with a residual scatter of < 2 subclasses. Both indices become 
independent of spectral type for spectral types earlier than -Ml 
and were used only for determinations between Ml and LO. 

Q and IH2O were calculated for the flux-conserved spectra 
of all targets. The two spectral type estimates, (SpTm)^ and 
(SpTm)/^20 which typically agree to within three subclasses, 
were averaged for each target. The subsequent detailed inspec- 
tion of these and another +three subclasses was usually sufficient 
to find the best match between the spectral templates and the tar- 
get component spectra. Targets with spectral types earlier than 
Ml were determined without this preselection. 

The final best-fitting template was selected by analyzing the 
absolute and relative strength of multiple Teff -dependent absorp- 
tion features in H and /T-band (e.g. Nal, Cal, All, Mgl, Si I, 
see Table [3]) by visual comparison of the continuum- subtracted 
target spectra with the selection of IRTF spectra. The result- 
ing spectral type uncertainties are typically between 0.5 and 2 
subclasses, estimated from the range of possible spectral type 
matches. Attempts to use an automatic spectral type determina- 
tion from line ratios did not return similarly good results. This is 
due to the weak strength of most absorption lines at the low spec- 
tral resolution ofR^ 1400 and the comparably uncertain contin- 
uum determination in the near-IR spectra. 

The template with the best spectral type was modified by 
adding artificial reddening and veiling in three iterations, each 



F\X) 



^phsph(^) 



+ k 



(4) 



(cf. iPrato et al]|2003l) . The best-fitting parameter combination 
and all un certainties were derive d through x^ minimization (for 
details see 'Daemgen et al.l l2012h . 

Owing to the relatively wide wavelength range of (-1.45- 
2.43 |LLm) of these observations, we tested whether the veiling 
can safely be assur ned to be independent o f A. As suggested by 
observations of e.g. lMuzeroUe et al] (l2003l) . veiling can take the 
shape of a blackbody at the dust sublimation temperature. To 
simulate this, the constant k was multiplied with a blackbody 
curve of temperature 1500 K and the fitting was repeated. The 
new best fits were at best indistinguishable from the fitting with 
constant excess, but mostly visibly worse, and the best x^ was 
also typically larger. Consequently, k was assumed to be constant 
for the derivation of vk- The resulting best-fitting spectral types 
A^y^^ and vk are listed in Table IH 

Best-fit templates typically fit the target spectrum very well 
over the wavelength range between ~1.7 and ~2.43 jam. The 
wavelength range between -1.5 and -1.7 jam is only fit well in 
-70% of the cases. Part of these deviations can be explained by 
uncertainties in the spectral extraction routine, which are larger 
at short wavelengths due to a lower Strehl ratio. The majority 
of deviations, however, are most likely due to broad water ab- 
sorption features at 1.38 jam and 1.87 |Ltm, which produce a char- 
acteri stic triangular co ntinuum shape in young late-type spectra 
(^e.g. JSato et al.|[l992l) that is not as pronounced in dwarf spec- 
tra because of their higher surface gravity. That surface grav- 
ity is indeed the reason for the deviant continuum at < 1.7 jam 
is demonstrated in Fig. O which shows that target components 
with the strongest bias toward the giant locus (i.e., with the low- 
est log g, intermediate between that of dwarfs and giants) are 
almost exclusively stars that strongly deviate from dwarf spectra 
at < 1.7 |Ltm. To estimate the impact of a poor fit at wavelengths 
< 1.7 |Ltm on the parameter extraction, test fits were run in a lim- 
ited wavelength range of 1.7-2.43 jam and compared with the 
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fitting over the whole range. The changes in the measured Ay 
and rK values were insignificant (|AAy| < 0.3, |Ar^| < 0.01) with 
larger uncertainties, however, than when derived from the whole 
available wavelength range. Accordingly, dwarf templates were 
fit to all target components in the whole 1.45-2.43 |Ltm range. 

The early spectral type (B9) of T 41 A is outside the range 
covered by the ^ectral library in use. This value was copied 
from the simbac^ database. No additional spectral values were 
derived for this component since discussing it is beyond the 
scope of this study. No spectral type could be derived for the 
brown dwarf candidates T 14 B and CHXR 73 B since the S/N of 
both spectra was too low to use spectral features for the deter- 
mination of the spect ral type. Assumi ng a spectral type of M9.5 
for CHXR 73 B (Luh manet aDl2006h . however, spectral fitting 
could be ran, resulting in an extinction value of A^^^ = 1.6 ± 0.7 
and veiling r^ = 0.02 ± 0.17. The spectrum of T 33 A shows 
signs of strong extinction and veiling. From the shape of the 
continuum, it could be inferred that its spectral type is earlier 
than or equal to M3. Due to the strong veiling, however, photo- 
spheric featur es were not stron g enough to further refine the SpT 
measurement. iLuhmanI (l2007h list T 33 A with a spectral type of 
G7, which the present data can neither confirm nor reject. The 
spectra of the components of ISO 126 and T 3 1 suff'er from high 
veiling or extinction, or from both. Spectral fitting for these stars 
was impossible due to their weak absorption features. 

We compared our results with spectral types and extinctions 
derived by Brandner & Zinnecker. (1997) . who observed the in- 
dividual components of T27, T31, and T39 with optical pho- 
tometry and spectroscopy. Their component spectral types are 
mostly identical to our NIR-derived values or agree to within 1 
(T 31 B) or 2 (T 39 A) subclasses. Their extinction values, how- 
ever, are lower by on average 0.6 mag for all target components 
except T39A. While this discrepancy can be explained by the 
subsequently assumed spectral type for this component, our sys- 
tematically higher values might be a result of the insufficiently 
known dust properties and thus extinction laws in regions of star 
formation and circumstellar disks, which makes an accurate ex- 
trapolation of vi sual extinctions fr om near-IR measurements dif- 
ficult (see, e.g., Nishiya ma et al .1 12009). The relative extinction 
measurements (Sect. 14.1.31) are not aff'ected by these systemat- 
ics. 



3.2.1 . Photometric vs. spectroscopic extinctions 

The spectroscopic measurements provide an independent esti- 
mate of the interstellar extinction in addition to the photomet- 
ric derivation in Sect. 13.11 Fig. [7] shows a comparison between 
both sets of extinction values. While the scatter around the line 
of equal extinctions is relatively high - deviations of 3 mag 
and more do appear - no systematic diff'erences could be no- 
ticed. There is also no systematic diff'erence between accreting 
components and non-accretors visible. The agreement between 
photometric and spectroscopic values indicates that the other 
spectroscopically inferred parameters (spectral type, veiling) are 
also on average compatible with the photometric results. Since 
the photometric determination of extinction, however, relies on 
shorter and thus more extinction- sensitive wavelengths (J-H 
color) than our spectroscopy (H-\-K), we used the photometri- 
cally derived values of all target components in the following. 




Ar^ [mag] 

Fig. 7. Spectroscopic versus photometric extinctions of all Cha I 
multiple star components for which both quantities could be 
measured. Red filled circles mark accreting components, all 
other components are marked with open circles. The dotted line 
shows equal extinctions. 



3.3. Effective temperatures and luminosities: the HR-diagram 

Spect ral types were convert ed to eff'ective tempe ratures accord- 
ing tolSchmidt-KaleJ (Il982h (earlier than MO) and lLuhman et al.1 
("2003") (later than or equal to MO). The luminosity of each tar- 
get component was estimated from its abs olute Mj magnitude 



|http : //simbad . u- strasbg . f r/simbad/| 



(assuming a distance to Cha I of 160+15 Dc:IWhittetet al.lll997h 
and bolometric corrections BQ (Hartigan et al .111 994 *). The lumi- 
nosity uncertainties were calculated from the propagation of the 
uncertainties of all input parameters; they are dominated by the 
uncertainty in the distance measurement. Luminosities, eff'ective 
temperatures, and stellar radii are listed in Table Hand are used 
to compile the Hertzsprung-Russell diagram shown in Fig. [8] 

3.3.1 . Stellar masses and ages 

The HR-diagram enables us to derive ages and masses of 
all components when comparing target locations with model 
isoch rones. To derive th e age and mass values in Table IH mod- 
els of ISiess et al.l (12000 ) were used. Because all derived values 
depend on the choice of mode l for the applied isochro nes, a com- 
parison of the same data with lPalla & Stable^ (Il999h isochrones 
is shown in Fig. [9l Qualitatively, there seems to be a slight 
bias toward younger ages in these latter models, because the 
sequence of low-c olor excess targ et components follows the 
3Myr isochrone of iPalla & Stable j closely, while the sam e tar- 
gets lie mostly below the 3 Myr isochrone of the Siess et al.l 
(2000) models in Fig. [8] (see Sect. 14.1.3] for a discussion of the 
derived component ages). Moreover, masses below -0.2 M© ap- 
pear on average to be slightly less massive when derived with the 
iPalla & Stabler' models instead of those by Siess et al. ("2000 1). In 
the following, the derived ages and masses from the Siess et al.l 
isochrones were u sed to guarantee comparability with the ONC 
results from .Daemgen et al.l (l2012l) . where the same model was 
applied. The differences between the inferred parameters from 
both sets of isochrones are probably smaller than their intrin- 
sic systematic uncertainties. While most results were only used 
for order of magnitude calculations, the inferred mass ratios, 
which is discussed in Sect. 14.2.21 have been shown to be quite 
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Table 4. Properties of individual binary and tertiary components 



o 






a 
X 



H 
H 



Oh 



& 



i 



Comp. 



SpT 



Teff 



[mag] 



rK 



[Z.o] 



[Myr] 



[Me] 



7?, 



[A] 



disk 
prob. 



log(Lacc/Lo) 



Mace 

[lO-'^Meyr^ 



B53 


A 


M2 + 1 


3560 + 150 


0.1+0.4 


0.00+0.09 


0.19+0.04 


4.o±t: 

8-0 ±3:2 


0.37 + 0.10 

90 + 0:07 
^■^^ - 0.05 


1.15+0.05 


<0.91 


0.00 


<-2.44 


<0.45 




B 


M4 + J;5 


3270 + 150 


0.2 + 0.4 


0.00+0.07 


0.05+0.01 


0.71+0.03 


<1.44 


0.00 


<-2.96 


<0.16 


CHXR 9C 


A 


Ml.5 + 1 


3630 + 150 


0.0 + 0.5 


0.04+0.16 


0.20 + 0.03 


4 4 + 5-1 

.3..3 + J0 


^■^' -0.07 


1.14 + 0.05 


<1.03 


0.00 


<-2.34 


<0.51 




BaBb 


M2±1.5 


3560 + 180 


0.1+0.5 


0.05+0.16 


0.22 + 0.03 


1.24 + 0.06 


<0.87 


0.00 


<-2.41 


<0.57 


CHXR 15 


A 


M5 + 1 


3130 + 140 


0.2 + 0.7 


0.00+0.15 


0.06 + 0.01 


4 5+2-5 


0.16 + 0.05 


0.87+0.04 


<1.01 


0.01 


<-2.88 


<0.29 




B 


M5 + I-5 


3130 + 170 


0.1+0.6 


0.00+0.13 


0.04 + 0.01 


0.14 + 0.06 


0.71+0.04 


<1.13 


0.03 


<-3.10 


<0.16 


CHXR 28 


AaAb 


K7 + ;-5 


4060+290 


2.7+0.4 


0.02+0.13 


1.2 + 0.3 


1.3+2-5 


0.75+0.25 


2.21+0.16 


<1.01 


0.00 


<-1.51 


<3.64 




B 


M2.5 + ;3 


3490 + 110 


0.4 + 0.6 


o.i3±S:i« 


0.42 + 0.09 


1.6 + 0.4 


0.36 + 0.04 


1.78+0.06 


<1.03 


0.01 


<-2.01 


<1.93 


CHXR 47 


A 


K5.5 + ;5 


4210+390 


5.8+0.7 


0.19 ±p 


0.72 + 0.14 


4.5 ±^-5 


0-95 ±S:S 

0.80 + J;20 


1.60 + 0.15 


<1.76 


0.17 


<-1.56 


<1.86 




B 


K7±;| 


4060+210 


5.2 + 0.8 


0.47+0.09 


5.5 ±|5 


1.39+0.07 


<1.78 


0.02 


<-1.81 


<1.08 


CHXR49NE 


A 


M2.5 + ; 
M3.5 + f 


3490 + 110 


0.0 + 0.4 


0.00+0.03 


0.20 + 0.04 


3-4 ±;| 


0.34 + 0.06 


1.22 + 0.04 


<1.07 


0.01 


<-2.46 


<0.50 




B 


3340 + 110 


0.0 + 0.3 


0.00+0.04 


0.15+0.03 


0.27+0.06 


1.15+0.04 


<1.12 


0.00 


<-2.45 


<0.60 


CHXR 68 


AaAb 


MO + 1 


3850 + 180 


0.3+0.5 


0.00+0.11 


0.44 + 0.10 


3.o±t-5 


46 + 0:10 


1.50 + 0.07 


<1.40 


0.04 


<-1.74 


<1.88 




B 


M1±J.5 


3710 + 110 


0.0 + 0.2 


0.00+0.05 


0.12 + 0.02 


12±j^^ 


0.83+0.02 


<0.82 


0.00 


<-2.79 


<0.12 


CHXR 71 


A 


M3±?-5 


3560+270 


2.3+0.6 


0.01+0.15 


0.11+0.03 


7.0 + 8-0 
13 +f 


0.37+0.05 


0.90 + 0.03 


<0.88 


0.00 


<-2.77 


<0.16 




B 


M5 + 1 


3130+270 


1.5+0.6 


0.04+0.18 


0.018+0.003 


013, 0.08 


0.44 + 0.04 


<1.05 


0.03 


<-3.51 


<0.04 


CHXR 73 


A 
B 


M2.5 + ;-5 


3490 + 180 


1.6 + 0.7 


0.02+0.17 


0.09+0.07 


8 + 50 


0.32 + 0.11 


0.83+0.04 


<0.98 

(<4.48) 


0.00 
0.03 


<-2.91 


<0.13 


CHXR 79 


A 
B 


M2 + 1 

M4 + f 


3560 + 150 
3270 + 110 


9.6 + 0.4 
6.4 + 1.0 


noo , 0.25 
^■^^±0.16 


0.36 + 0.07 
0.035+0.007 


2-0 ±Q? 
12 ±55 




1.57+0.07 
0.59+0.02 


0.92 + 0.21 
<1.96 


0.74 
0.21 


-1.53+0.24 
<-2.98 


4.88±3-«3 
<0.13 


Hn4 


A 


M4 + ;-5 


3270 + 110 


3.6 + 0.5 


0.01+0.12 


0.10 + 0.02 


4.2 + 1.0 


22 + 0-0^ 

^■^^ -0 05 


1.00 + 0.03 


<1.61 


0.00 


<-2.42 


<0.69 




B 


M3 + 1 


3420 + 150 


3.5+0.5 


0.00+0.10 


0.11+0.02 


CO ,5.0 
D.Z + Q7 


0.30 + 0.08 


0.95+0.04 


<1.72 


0.00 


<-2.39 


<0.52 


Hn21 


A 


M3.5 + ;-5 
M6.5 + ^5 


3340 + 110 


2.6 + 0.7 


0.02+0.14 


0.11+0.02 


4 5+2-0 


0.26 + 0.05 


0.99+0.03 


<0.99 


0.18 


<-2.81 


<0.24 




B 


2940+200 


0.0 + 0.6 


0.00+0.14 


0.030+0.006 


o.o8±S:;o^ 


0.66 + 0.04 


<1.05 


0.01 


<-3.22 


<0.20 


ISO 126 


A 


M0.5 + i'5 


3780+400 


>4 


»1 


0.50 + 0.09 


2.0 + ?! 


0.52 ±s-s 


0.47+0.05 


(<0.75) 


0.00 








B 


MO.5 + 1-5 


3780+400 


>4 


»1 


0.090+0.006 


(23±T8) 


0.50 + 0.2 


0.70 + 0.07 


(<0.65) 


0.00 






T3 


A 


M1.5 + J-J 
M3 + i^ 


3630+320 


2.2 + 1.1 


J 3 + 1.8 
0.18+§| 


0.50 + 0.09 


l-6±?-^ 


O.42 + S-51 


1.79+0.16 


6.52 + 0.28 


1.00 


-0.62+0.05 


40.9+51-4 




B 


3420+250 


0.5+0.7 


0.10 + 0.02 


5.8 ±g 


00 +0.13 
^■^^±0.10 


0.91+0.07 


<1.37 


0.00 


<-2.45 


<0.45 


T6 


A 


Kl + 3 


5080+490 


0.0 + 0.9 


1.5 +J1 


1.6 + 0.4 


10 + 12 

4-2 ±i§ 


1.4 + 0.2 


1.65+0.16 


<0.85 


0.00 


<-1.37 


<2.01 




B 


M6 + J5 


2990+320 


4.6 + 1.2 


0.54+0.14 


0.06 + 0.01 


10 + 0-11 
U.iU + 07 


0.89+0.09 


<2.00 


0.00 


<-2.55 


<1.00 


T14 


A 


K7±;:3 


4060 + 180 


0.1+0.5 


0-33 ±S- 


0.67+0.15 


3.0 + 3-0 


0.78+0.20 


1.66 + 0.07 


5.86 + 0.33 


0.00 
0.00 
1.00 


-0.90+0.04 


10.7+3.0 


T26 


A 


GO + 1-5 


6030+290 


6.6 + 0.9 


0.18 ±SiJ 


4.7 + 1.0 


7 1+3-0 


1.7+0.2 


2.00 + 0.10 


5.62 + 0.30 


-0.50+0.15 


14.9+6.4 




BaBb 


M5 + 1-0 


3130+290 


0.9+0.8 


0.00+0.16 


0.11+0.02 


0-18 ±SS 


1.13+0.11 


<1.18 


0.01 


<-2.49 


<0.81 


T27 


A 


Ml.5 + 23 


3630 + 180 


2.9+0.6 


0.10+0.19 


0.14 + 0.03 


-:! 


o-4o±S:;j 


0.96 + 0.05 


1.83+0.22 


1.00 


-2.07+0.25 


QO + 0.65 
^■°^-0.70 




B 


M3 + 1 


3420 + 150 


2.5+0.7 


0.00+0.12 


0.09+0.02 


0.29+0.07 


0.87+0.04 


<1.00 


0.01 


<-2.85 


<0.17 


T31 


A 


K7±?5 


4060+460 




(»1) 


1.6 + 0.3 


1.0±3.° 




2.55+0.29 


(9.35+0.19) 


1.00 


>-1.51 


>4.20 




BaBb 


MO + 1.5 


3850+290 




(~0) 


0.76 + 0.10 


1-5 ±5l 


1.97+0.15 


(<0.87) 


0.01 






T33 


A 


<M3 






(»1) 










(<0.84) 


0.00 








BaBb 


M0.5 + [5 


3780+210 


1.5+0.5 


0.00+0.12 


0.57+0.13 


l.i±lj 


52+0-18 
U.DZ + QJ4 


1.77+0.10 


<0.99 


0.00 


<-1.83 


<2.01 


T39 


A 


M2.5+0.5 


3490+70 


0.2 + 0.4 


0.00+0.08 


0.25+0.05 


2-7 ±2:1 

4.0 ±; 

4.1 ±\l 


0.36 + 0.05 


1.37+0.03 


<1.02 


0.00 


<-2.17 


<1.03 




B 


M2 + 0.5 


3560+70 


0.2 + 0.4 


0.00+0.07 


0.20 + 0.04 


0.38+0.05 


1.17+0.02 


<1.05 


0.00 


<-2.33 


<0.58 




c 


M2.5+0.5 


3490+70 


0.0 + 0.3 


0.00+0.07 


0.14 + 0.03 


0.34 + 0.05 


1.04 + 0.02 


<1.12 


0.01 


<-2.42 


<0.47 


T41 


A 


B9^ 


10500 














<0.36 


0.00 








B 


M3.5+0.5 


3340+70 


0.0 + 0.2 


0.00+0.07 


0.11+0.08 


4 5 + " 


0.26 + 0.07 


0.99+0.02 


<1.28 


0.00 


<-2.44 


<0.55 


T43 


A 


M3.5 + 1-'' 


3340+290 


3.4 + 0.7 


0.44+0.16 


0.29+0.06 


4.0 ±« 


0^08 + S|^ 


1.62 + 0.13 


<1.07 


0.00 


<-2.34 


<0.99 




B 


M7 + 1-0 


2880+360 


4.5+0.8 


0.00+0.16 


0.05+0.01 


0.91+0.11 


<1.07 


0.00 


<-3.00 


<0.46 


T45 


A 


Ml + {'3 


3710+200 


2.8 + 1.2 


1.5+2-2 
22 + 0-19 


0.81+0.15 


1 1+0.6 

33+1^ 


0.49+0.12 


2.18+0.12 


12.62+0.47 


1.00 


-0.26+0.10 


98+35 




B 


M5±i, 


3130+210 


1.6 + 0.9 


0.10 + 0.02 


0-18 ±S:S 


1.08+0.07 


<1.28 


0.01 


<-2.62 


<0.58 


T51 


A 


K7 + 3 


4060+590 


1.1+0.7 


0-48 +S-S 


0.76 + 0.17 


2.5 ±{3, 


0.78+0.40 


1.77+0.26 


<1.08 


0.01 


<-1.55 


<2.56 




B 


M2.5 + 1.5 


3490+220 


3.8+0.8 


0.10+0.20 


0.10 + 0.02 


7.2 ±52 


37 +0-12 


0.85+0.05 


1.11+0.38 


0.68 


-2.67+0.47 


0.23+0.45 


T54 


A 


G8 + 3 


5507+440 


0.2 + 0.5 


0.00+0.10 


2.0 + 0.4 


16±30 
9.O + J63 


1.3+0.2 


1.52 + 0.12 


<0.87 


0.48 


<-1.69 


<0.96 




B 


K7 + 1.5 


4060+310 


0.0 + 0.5 


0.04+0.13 


0.35+0.08 


0.80 + Jf 


1.20 + 0.09 


<1.07 


0.00 


<-2.23 


<0.35 



Notes. *^^^ Spectral type from Simbad. ^^^ Values in parantheses mark lower limits since no veiling could be measured. ^^^ Masses <0.1 M© lie outside the scope of the lSiess et alJ J200Qh models. 
Values are estimated through extrapolation. 
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5.4 



Fig. 8. HR diagram with the isochrones of ' Siess et alJ (l200d ). 
Filled red circles show spatially unresolved components, ei- 
ther close binary components of visual triple stars or spec- 
troscopic binaries. Luminosities of these targets probably ap- 
pear to be higher than single stars at the same Tqq. Open blue 
squares show targets with strong near-IR excess (E(Ks-L') > 0.5 
or H-K^ > 0.6) and filled blue squares all other components of 
the ChamaeleonI sample. In the lower left corner we present the 
average random error and the luminosity uncertainty introduced 
by 0.2 mag variability in /-band (estimated from measurements, 
see Sect l2.2.2T) . 




G5 KO 

Spectral Type 



Fig. 10. Equivalent widths of Bry absorption lines as a func- 
tion of spectral type measured from the IRTF dwarf template 
sequence. The continuous red line shows our best fit (Eq.O. 



ity ( 'Fischer et al.ll201 ll) . Additionally, extinction might be over- 
estimated, as discussed in Sect. 14.1.31 which can also lead to 
higher derived luminosities. These eff'ects explain why the ten 
stars with the youngest derived ages are all either unresolved bi- 
naries or stars with strong color excess. 
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3.5 



3.4 



Fig. 9. Same as Fig. [H but with isochrones of |Palla& Stable j 
([l999h . 



insensitive to the choi ce of isochrone (e.g., ICorreia et al.ll2013l: 
iLafreniere erani2008h . 

As expected. Fig. [8] shows that unresolved binary stars have 
higher-than-average luminosities compared with individually 
observed components at the same Tq^. Depending on the bright- 
ness ratio of these close components, the derived ages can ap- 
pear artificially young. The same is true for stars with strong 
near-IR excess. While excess is typically assumed to be weak 
in the /-band (see, e.g., Meyer etal. 1997), stars with overall 
strong excess emission can exhibit enough excess flux in J- 
band to significantly increase the estimated bolometric luminos- 



3.4. Accretion 

Accretion activity was inferred from the existence of Bry 
(21 665 A) emission. To assess the strength of this line, we used 
the high-S/N spectral extraction. 



3.4.1. Equivalent widths ^Bry 

Equivalent widths were integrated in a 5 6 A- wide interval around 
the centered Bry peafl The continuum and noise background 
were estimated from regions shortward and longward of the 
line's central wavelength. Bry measurements were performed 
on the calibrated spectra, correcting for veiling by applying 
W^ = (l-\- vk) X Ty^easured (Daemgen^etaD[l012). To correct for 
the intrinsic Bry absorption of spectral types earlier than Ml, we 
measured the equivalent width in the same wavelength range for 
the spectral sequence of dwarf stars in the IRTF database. The 
correlation between W^ry and spectral type is shown in Fig. [TOl 
and was fit by a quadratic polynomial 



^photosph^3p^^ = 



(0.16 ± 0.05) - (0.16 ± 0.01) X SpT 
+ (0.0030 ± 0.0004) X SpT^ , 



(5) 



valid in a spectral type range between FO and Ml, with spectral 
types (SpT) measured as M0=0, M5=5, K0=-8, and so fortlfl. 
The absorption for spectral types later than Ml is low and no 
correction is applied here. Although a rather large correction is 
probably necessary, spectral types earlier than FO remained un- 
corrected because the spectral template sequence does not ex- 



^ Absorption features are defined to have negative equivalent width 
values and emission is positive. 
^ The adopted classification omits spectral types K8 and K9 
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tend to the earliest spectral type of this sample (T41 A, B9). The 
correction terms, tabulated in Table |5] were applied as 



w^^ = w. 



w 



photosph 
Bry 



(6) 



resulting in equivalent widths W^^^ that are free of veiling and 
photospheric absorption. We assumed that the measured equiv- 
alent widths are a good indicator for the accretion, because Bry 
has been found by lAntoniucci et al] (1201 ll) to be well correlated 
with the accretion luminosity of young low-mass stars, better 
than the other examined tracers Ha, Call, [O I], and Pap. When 
the uncertainty of the equivalent width measurement was larger 
than the derived accretion signal, upper limits were derived from 
a noise measurement in the continuum surrounding Abyy- The 
resulting fully corrected equivalent widths or upper limits are 
Hsted in Tabled 

To infer the presence of an accretion disk around all compo- 
nents, the disk probability (Daemgen et al. 2012) was calculated 
for each target component from the measured W^^^, its uncer- 
tainty, and the local noise. The values are listed in Table IH 



3.4.2. Line luminosities and mass accretion rates: Lsry, ^acc, 

&Macc 

Line luminosities were derived from the flux-calibrated and 
dereddened spectra and were converted to accretion lumi- 
nosities using the emprical correlatio n log(Lacc) = ( }-^^ ^ 
0.19)log(LBrr/Lo) + (4.43 ± 0.79) from lMuzerolle et al.l (Il998h . 
Mass accretion rates Mace were calculated for all targets with 
significant accretion according to 



Mace = 



^acc ^* 



Ri, 



GM^ \ 7?in — R^ 



(7) 



(iGullbring e t al.|[T998 ^), assuming /?in ~ 5/?*, with stellar radius 
R^ derived from Stefan's law using the individual stellar lumi- 
nosities and eff'ective temperatures. Upper limits for target com- 
ponents with no detected Bry emission were determined from 
the minimum-emission feature detectable with 3cr significance 
over the noise in the vicinity of the line. All inferred star and 
disk parameters are listed in Table |4] 

3.5. New companion candidate to T33B 



T 33, also known as Gla ss I (lGlassiri 979), was discovered to have 
a binary companion by IZinnecker et al. (1988). Our new obser- 
vations show for the first time clear indications tha t the sec- 
ondary component T 33 B, which was reobserved by Ghe z et all 
(1997) with a separation of 2r5±a'5 and PA=(284±5)°, it- 
self consists of two stellar components (in the following called 
T 33 Ba and T 33 Bb; see Fig. O. To test whether the newly de- 
tected tertiary component is a real companion, we considered 
component colors, the surrounding stellar density, and proper 
motion of the components. 

The dereddened colors of t he two compjpnents of T33B 
agree roughly with PMS colors ("Luhman et al.l l201Ql) of K7 and 
Ml spectral type. Alternatively, a spectral type of K4 can be 
assigned to the two c omponents assuming that they are giants 
(iBesseU & Brettlll988 ). The spectral type of the combined fight 
was derived to be MO.S^j^. Because the total light from the 
T 33 B system is dominated by its brighter component, this con- 
firms the PMS nature of the primary of the T33B subsystem, 
since a spectral type as early as K4 can be excluded by the de- 
tected photospheric features in H-\-K. 
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Fig. 11. Position of T 33 Ba (right group of points) and T 33 Bb 
(left group) in RA and DEC relative to T 33 A at four epochs. 
The position of T 33 Ba & Bb at the earliest epoch has a com- 
parably large uncertainty since both components are spatially 
barely resolved. 



iLafreniere et al.l (l2008l) evaluated the probability of chance 
alignments of multiple components based on the density of stars 
in the 2MASS catalog in the closest 15' around each star of 
their sample. For all but the widest components they calculated a 
probability of < 10""^ for chance alignments. Since this includes 
the 2! 4 wide T 33 A-hB system, the new component located only 
O'.'l from T 33 B must have a considerably lower probability for 
a chance alignment, making a physical association very likely. 

To fully ascertain whether Bb is a background object or 
physically bound to the T 33 system we used proper motion mea- 
surements of the individual components. Two additional epochs 
of NACO observations taken during the past six years were avail- 
able in the ESO archivCl both of which resolve T 33 B into two 
components. They use a similar instrumental setup as the current 
observations (^s filter, S13 camera). Together with the proper 
motion of T 33 A o f ^r a = -40 mas/yr and ^dec = +33mas/yr 
(iTeixeira et al.ll200Q) , the common proper motion is clearly con- 
firmed: within the -five years between the first and last of the 
observed epochs, both T 33 Ba and T 33 Bb move less than O'.'OS 
in RA and DEC and in diff'erent directions relative to T33A 
(Fig. [IH). The expected motion of T 33 A of ^ra=0'.'20 and 
^DEC=0'-'17 in the same time span renders the two components 
comoving with T 33 A, and thus T 33 is very likely a physically 
bound triple stellar system. 

The residual motion between the components T33Ba and 
Bb is probably orbital motion. An estimate of the shortest or- 
bital period using Newton's and Kepler's laws under the as- 
sumption of a circular orbit and that the last observation was 
at maximum elongation returns Pniin=102±15 yrs. Classification 
of T 33 BaBb as a binary that has been hiding its binary nature 
through aligning its components in an inclined orbit during the 
past few years explains why its binary natur e has not been dis- 
covered before and why lNguyen et al.l(l2012b observed a deviant 
radial velocity of T 33 B from the rest of the cluster while no sig- 
nal of spectroscopic binarity was detected. 



8 2006-03-25 [076.C-0579(A)], 
ihttp : //archive . eso . org/| 



2008-02-20 [080.C-0424(A)], 
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Table 5. Equivalent widths of Bry in dwarfs 



SpT 



GO 
Gl 
G2 
G3 
G4 



Wl 



rphotosph 



[A] 



-4.00+0.23 
-3.74+0.21 
-3.48+0.19 
-3.23+0.18 
-2.98+0.17 



Wl 



photo sph 



SpT 



Bry 
[A] 



Wl 



photo sph 



SpT 



Bry 
[A] 



G5 
G6 
G7 
G8 
G9 



-2.74+0.15 
-2.51+0.14 
-2.28+0.13 
-2.05+0.12 
-1.84+0.10 



KO 
Kl 
K2 
K3 

K4 



-1.63+0.09 
-1.42+0.09 
-1.22+0.08 
-1.03+0.07 
-0.85+0.06 



Wl 



SpT 



•photo sph 
Bry 
[A] 



K5 
K6 

K7 
MO 
Ml 



-0.67+0.06 
-0.49+0.05 
-0.32+0.05 
-0.16+0.05 
0.00+0.05 



4. Discussion 

4.1. Stellar parameters and sample biases 

4.1 .1 . Full, extended, and restricted samples 

The full sample of 26 multiple (= 52 individually detected) ob- 
jects comprises targets from a broad range of binary separations, 
component spectral types, and observational quality. To ensure 
the comparability and validity of the results, we defined in addi- 
tion to the full sample of components with individually obtained 
spectroscopy (irrespective of the fact that they might be close bi- 
naries) two new samples, the extended sample and the restricted 
sample. 

The full sample consists of all components "as observed". 
This disregards that some of these components may be unre- 
solved visual or spectroscopic binaries. The extended sample, 
in contrast, is composed of all known individual components. 
Spectroscopic binaries as well as known visual components that 
remain unresolved by the present study each enter the sample 
as two individual components. Since not all individual stellar 
parameters are known, this sample will only be used for num- 
ber statistics such as the disk frequency (Sect. 14.2.11 ), where the 
presence of disks in unresolved binaries can be inferred from 
the unresolved data. This sample grants the most pristine view 
to the data, since no additional restrictions to the original com- 
ponent sample are applied. The extended sample consists of 62 
components (19 binaries, 4 triples, and 3 quadruple stars). 

The restricted sample was composed to eliminate observa- 
tional and compositional biases by applying the following selec- 
tion criteria, (a) Spectral types must be in the range for T Tauri 
stars, i.e., later than GO. (b) All target components must be stel- 
lar. It is not clear whether the formation process of brown dwarfs 
is identical to that of stars. To exclude the eff'ects of possibly dif- 
ferent initial conditions, all brown dwarf candidates in the sam- 
ple were excluded, (c) The information must be attributable to 
one stellar component. A few stars in the sample could not be 
resolved with our adaptive optics spectroscopy or are suspected 
spectroscopic binary candidates. This requirement limits the re- 
stricted sample to binaries with projected separations of > O'.'IS 
or > 25 au at the distance of Cha I. While these restrictions mean 
that some of the binary components are "singled", i.e., only one 
component enters the sample, it reduces the slight preference 
for primaries to have earlier spectral types than secondaries (see 
next paragraph and Fig. [121) since early-type primaries and late- 
type secondaries are removed. 

The restricted sample comprises 43 components. If binary 
parameters are to be investigated with respect to the restricted 
sample, only binaries with both components in the restricted 
sample are considered. This restricted binary sample consists 
of 17 binaries. 
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Fig. 12. Histogram of spectral types divided into primary stars 
(solid orange) and secondary stars (black outline). The spec- 
tral types of secondaries are on average 1.5 subclasses later. 
Unresolved close binaries in triple systems are plotted as one 
component. 



4.1.2. Spectral types 

As can be seen in the spectral type distribution in Fig. [121 the 
median spectral type of primary stars (Ml. 5) is earlier than that 
of the secondaries (M3). This is true for the full and restricted 
samples. To test whether the subsample of binary stars is rep- 
resentative for the population of Cha I members, we calculated 
the probabilities that the spectral types of primaries and secon- 
daries are dra wn from a com mon parent distribution with single 
stars of Cha I (lLuhmanl2QQ4 . 2007), respectively. We limited the 
single-star sample to V<16.7 mag to account for the fact that the 
primaries were selected to serve as SINFONI adaptive optics ref- 
erence stars. The resulting KS probability of 92% indicates that 
the distribution of primary star spectral types in this sample can- 
not be distinguished from the single star distribution. Secondary 
stars, however, are not subject to this selection bias and agree 
with the unconstrained sample of single stars (KS probability 
19%). 

The apparent diff'erence between the primary and secondary 
spectral type distributions in Fig. [12] by definition reflects the 
mass difl'erences between the sets of primary and secondary 
components. This can produce biases when comparing the pa- 
rameters of both populations since many parameters, such as 
stellar evolutionary timescales, are mass-dependent with the 
more massive components typically evolving faster than the less 
massive stars. The resulting difl'erences of the inferred parame- 
ters are discussed in the subsequent sections when they are sus- 
pected to occur. 
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4 6 

Av(primary) [mag] 

Fig. 13. Primary vs. secondary extinctions Ay ^\ Target bina- 
ries with at least one component with strong near-IR color ex- 
cess E(K^-L') > 0.5 are marked by a red circle. The dotted line 
shows equal extinctions in both compononents. The three ma- 
genta stars (CHXR73, T26, T33) have at least one component 
without measured E(K^-U). From their very red color of H-K 
>0.6mag a strong near-IR excess can also be inferred for these 
binaries. 
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Fig. 14. Primary age vs. secondary age in logarithmic units. 
Open symbols are target binaries with unreliable age measure- 
ments because either their primary (triangles right) or secondary 
(triangles up) or both components (diamonds) are spectroscopic 
binaries or exhibit a strong color-excess E(K^-L')>0.5 (or 
H-K^ >0.6, if K^-U could not be measured). Filled symbols 
with error bars are binaries with spatially separated components 
with low color excess. The dotted line is the reference for coeval 
ages. 



4.1 .3. Relative extinctions and ages 

In Fig. [13] we compare the photometric extinctions in both com- 
ponents of each binary. While most of the extinctions are iden- 
tical for the components of the same binary, the multiples with 
significant diff'erences in component extinction are mostly host 
to at least one very red component with near-infrared excess 
E(K^-L') > 0.5. If extinction is dominated by the degree of em- 
beddedness of a binary in the molecular cloud, it is expected 
to be similar for both components of a binary. The presence of 
K^-L color excess apparently leads to a diff'erent photometric 
extinction value for the respective component. Since the spec- 
troscopically determined extinctions (see Table lU show no such 
correlation (the values scatter almost uniformly around equal 
primary and secondary extinctions, though with a larger spread 
due to the larger uncertainties of A^^^), additional systematic 
uncertainty of the photometric extinctions as a function of color 
excess may be present that do not aff'ect the spectroscopically de- 
rived extinctions. A possible origin can either be an overestima- 
tion of the slope of the extinction vector in the (J-H)-(H-K^) 
color-color diagram (Fig.|3]l or an underestimation of the slope of 
the CTTS locus (which can be caused, for example, by assuming 
filter system transformations that do not accurately describe the 
filters in use) the result would be a mostly positive correlation of 
Ay with the color-excess, as observed. In the following, we used 
the derived extinctions as measured and caution that a possible 
systematic uncertainty may be present with a median amplitude 
of 2 mag for high E(K^-U) components. If true, the luminosity 
of the respective components would be overestimated by on av- 
erage 0.21 dex, comparable in size to the random uncertainty of 
the luminosity values, resulting in an age estimate that is younger 
than those derived in Table|4]by A log r ^ 0.3 at M = 0.4 M©. 

Taking these eff'ects into consideration, the measured ages 
show that binary components are coeval. Fig. [14] compares 
the logarithmic isochronal primary and secondary ages from 
Table [4] Binaries with at least one component with strong color 



excess or an unresolved binary component are mostly located off" 
the sequence of equal-ages. This is due to overestimated lumi- 
nosities of components with strong continuum excess or over- 
estimated extinction and accordingly under-estimated ages. In 
contrast, all binaries with low-excess components are consis- 
tent with equal age components to within Icr. To test whether 
the observations are compatible with the simultaneous forma- 
tion of the two components of a binary, we compared the log- 
arithmic age diff'erences |AlogT| = llogTprim - logTsed follow- 
ing Kraus & Hillenbrand (2009). The rms scatter in |AlogT| is 
0.03 dex, which is compared with 10000 randomly paired sam- 
ples of the same size as the target sample, drawn from the pri- 
maries and secondaries with low veiling. Only 2% of the random 
samples have an rms scatter smaller than 0.03 dex. The binary 
sample is thus more coeval than the total sample of binary com- 
ponents with a significance of 2 3cr. 

The median inferred age of the subset of low-r^ compo- 
nents (excluding unresolved binaries) is 5.5 Myr, older than the 
2-3 Myr derived by Luhman (2007) for a large sample of Cha I 
members with models of Baraff'e et al. (1998) and Chabrier et al.l 
(20 0Qb. Even th e systematically younger ages when using the 
iPalla & Stahlei] ([1999) models (see Fig. [9]) produce a median age 
(4 Myr) that is older than the single- star value. Since it is un- 
likely that binaries are intrinsically older than single stars of the 
same region, this discrepancy must be due to either sample bi- 
ases or a systematic diff'erence in the derivations of luminosities 
and eff'ective temperatures. While part of the diff'erence might 
arise from the isochrones used, we found that the luminosities 
derived by iLuhmanI (l2007l) are systematically higher than those 
derived in the binary s ample. Pa rt of this diff'erence is due to 
the larger distance that iLuhmanl assumes. The increase of the 
distance from 160 to 168 pc requires a luminosity increase of 
A log L = 0.05 for all targets to match the observed /-band mag- 
nitudes. Furthermore, while a correctio n factor for unresolved 
binarity was applied by iLuhmanI (l2007h . he did not account for 
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possible excess emission in /-band, which can lead to compara- 
bly high luminosities, as we see for the high- veiling components 
in the HR-diagram (Fig. [8]). This can lead to an underestimation 
of the median age. In fact, if one includes all binary components 
of the present sample irrespective of their veiling into the age es- 
timatio n, a median age of 4.4 My r is derived froni thelSiess et al.l 
(12000') models and 2.6 Myr from'Palla &Stahler| (Il9m 

A thorough absolute age determination should be repeated 
with a larger sample than the one presented here. Nevertheless, 
the age determination serves to confirm the similar relative 
ages of the binary components as discussed above. Coevality 
and equal component extinctions a re in line with previous find- 
ings of binaries in Ch al (Brandner & Zinnecker 1997) and other 
star- forming regions (IWhite & Ghezl200ltlKraus & Hillenbrand! 
|2009|). 



4.2. The evolution of disks In Cha / binary components 

4.2.1. Disk frequency 

From the number of target components with K^-L' color excess 
that were used as an indicator for the presence of hot dust in 
the inner disk and those with significant Bry emission, which 
indicate ongoing accretion activity, we determined the frequency 
of dust and accretion disks around the components of binaries 
in Cha I and compared it with disks around single stars of this 
region. 
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Fig. 15. Near-infrared excess E{K^-L) for the primary and sec- 
ondary components of each binary with individual component 
photometry. Red filled circles and blue filled diamonds are bi- 
naries with accreting primaries and secondaries. Targets to the 
right and above the dotted lines have significant near-IR excess, 
indicating the presence of inner disk material. A typical error 
bar, dominated by the systematic uncertainty of -0.2 mag, is in- 
dicated in the upper left. 



Accretion disks. The fraction of target components harboring 
accretion disks was determined from the indivi dual disk proba- 
bilities in Table |4] as described in the appendix ofl Paemgen et al.l 
(l2012h . Some components have no good estimates for the veil- 
ing level. A correction of W^ry is not possible for these targets, 
and the accretion strength might be significantly larger than mea- 
sured. One of the reasons for an inaccessible veiling value is high 
veiling itself, since the line features used for the veiling estima- 
tion can be reduced in depth until they blend into the continuum 
noise. Excluding these targets may accordingly bias the sample 
against accretors. To avoid this bias, these targets were included 
with a neutral disk probability of 0.5. As a test, we also derived 
disk fractions assuming the limiting cases of 0% and 100% disk 
probability and found no significantly diff'erent results. In the 
following, we consider the extended and the restricted samples. 

The extended sample of all components returns an accretion 
disk fraction of F^^^ = 18^^%. This calculation is possible since 
all unresolved binaries except T 3 1 A show no Bry emission and 
can be treated as consisting of two non-accreting components. 
T 31 A is a suspected SB 2 and shows signs of accretion (Wery 
> 9.4 A). The uncertainty of being host to one or two accreting 
components is included in the calculation of the disk frequency. 

The restricted sample returns a very similar frequency of 
^restr _ 16+9% Bq^]^ Calculations include only the informa- 
tion from the Bry emission measurements. It is known that 
Bry returns a lower fraction of accretors than the H a feature, 
which, however, defines the classes of classical and weak-line 
TTauri stars. This is due to the ~100 times lower transition 
probability of Bry excitation compared to Ha B With the help 
of th e correction factor / = 0.125 derived in iDaemgen et al.l 
(1201 2) from the fraction of accretors detected in Ha and Bry 
and a number of A/^ = 43 components in the restricted sample, 
the number of CTTS without Bry emission is estimated to be 
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f/(l-f)x(F-N)^l. This increases the restricted accretion 
disk fraction to Flf^' = \9^_l%. 

This fraction is significantly lower than the accretion disk 
fraction of single s tars in Chal of F^'""^^^ = 44 ± 8% 
(iMohanty et al.|[2QQ5h . Assuming that binaries and single stars 
are bom at the same time in the star formation history of a clus- 
ter, this indicates that binaries with separations of -25-1000 au 
finish accreting on a shorter timescale than single stars. To in- 
vestigate this trend in more detail, we divided the sample into 
binaries with separations closer and wider than 100 au, close 
to the median separation of 119au of the sample. The result- 
ing fractions of F'^f^^^^ = 10^^^% and F'^f^^^^ = 1^}^% show 
that the accretion disk fraction is particularly low in the clos- 
est binaries of the sample, less than half of the wider sample. 
This corroborates theoretical expectations that the timescale for 
disk evolution is significantly shortened in binaries with sepa- 
rations <100au, if one assumes a correlat ion of disk lifetini e 
with the size of the (truncated) disk (e.g.. iMonin et al.ir2007h . 
The fact that the reduction of disk frequency in wide bina- 
ries is less pronounced but still significant indicates that disk 
lifetimes can be shortened even if the truncation of the outer 
disk is only weak. Th is agrees with the results from the ONC 
(Daemge n et al.l2012h . where the sample of binaries does show a 
reduction of accretion disk frequency although only separations 
of -100-800 au and no closer pairs could be examined owing to 
the comparatively large distance of the ONC. 

Inner dust disks. The frequency of inner dust disks g is mea- 
sured from the number of targets with significant K^-L excess 
in the nearly simultaneous photometry. A diagram of primary 
and secondary E{K^-L) is shown in Fig. [T2l As derived in 
Sect. 13.11 the systematic uncertainty of the K^-L excess mea- 
surements is about 0.2 mag. Accordingly, K^-L! emission is in 
significant excess over the pure stellar color for components 
with E(K^-U) > 0.2. There are 21^^^ components out of 55 
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with measured K^-L colors in the extended sample fulfill this 
criterion, and 20^^^ components in the restricted sample of 40 
components with measured colors. Uncertainties of the number 
counts were determined from targets for which E{K^-L) = 0.2 
is within the uncertainty limits. The resulting inner dust disk 
fractions are 3^^^ = 38+^^ % and W' 



50^^^ %. 



The values agree well with p revious measurements of the 
disk frequency in binaries of Cha I. iDamjanov et al.l (l2007h mea- 
sured the frequency of targets with 8 or 24 jam excess in 15 spa- 
tially unresolved binaries and 2 triples in Cha I and derived a 
nominal frequency of multiples with disks of g = 35^|3%. Their 
value, obtained from 10-300 au binaries, agrees well with the 
fraction of components with K^-L excess obtained from the ex- 
tended sample. And indeed, wh en applying the s election criteria 
of the restricted sample to the iDamjanov et al.l binary sample, 
the inferred hot dust disk fraction increases to 42%, since the 
fraction of binaries with separations < 25 au that show no signs 
of dust excess is high (four out of five). This is compatible with 
the fraction found for the restricted sample of binaries. 

The extended and in particular the restricted inner disk frac- 
tions of Cha I binaries are compatible with the average Cha I 
single star value of -47-55% (with weak fluctuations depend- 
ing on stellar mass) found from photometric excess in the 
Spitzer/IRAC bands between 3.6 jam and 24 |Ltm (iLuhman et al.l 
1200 8: Damjano v et al.l l2007l) . A significant reduction of disk 
presence compared to single stars is mostly visible in close bi- 
naries: out of 16 components that are part of spectroscopic or 
visual binaries with separations < 25 au, only 3 or 4, depending 
on whether the spectroscopic binary T 3 1 A has one or two disk- 
bearing components, were found with significant K^-L excess, 
a fraction of 19-25%. The inclusion of close binaries explains 
why the extended sample shows a lower disk fraction than the 
restricted sample. 



The correlation of inner disk presence and accretion. The 
fraction of inner dust disks around binary components is signif- 
icantly higher than the accretion disk fraction. Figure [16] com- 
pares measured Bry equivalent widths with the near-IR color 
excess. While all accreting components show significant color 
excess, a nominal fraction of 65% of the binary components with 
near-infrared excess shows no signs of accretion in Bry. This re- 
duces to 55% when accounting for the fact that ISO 126, the bi- 
nary with the two components with the strongest near-IR excess, 
has been found to be accreting in Ha studies ( Saff'e et al. 2003). 
This fraction is significantly higher than for Cha I single stars ac- 
cording to the study by .Damjanov et al.l (I2QQ7.) . who found 6 of 
21 stars with inner dust disks (inferred from [4.5] -[8.0] excess) 
without signs of accretion (Ha). 

This discrepancy of the accretion and dust disk presence has 
been observed not only for Cha I, but in several star- forming re- 
gions, making a stochastic explanation accretion variability un- 
likely. We conclude that the approximately exponential decay of 
disk frequency is more rapid for accretion disks (time constant 
Tacc ^ 2.3 M yr) than for hot inner dust disks (rdust ~ 3 Myr; 
iFedele et ani2010,) . Binary companions in Cha I apparently am- 
plify the effect that accretion disks disappear before dust disks 
around the same stars. Possible reasons are an inside-out disk 
evolution, UV irradiation from the star, or planet formation and 
migration. A detailed discussion will follow in a forthcoming 
paper (Daemgen et al. 2013). 
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Fig. 16. Bry equivalent width as a function of near-IR color ex- 
cess E(K^-L'). Filled symbols show accretion detections, open 
symbols are upper limits of Wbvy • Primaries and secondaries are 
denoted by circles and diamonds. Stars to the right of the dot- 
ted line have significant K^-L' excess. Targets without veiling 
measurement can have larger W^ry than shown here since the 
correction factor (1 + r^) could not be applied; they are marked 
with arrows. 
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Fig. 17. Accretion status as a function of binary separation. 
The gray shaded bars show the restricted sample, the combined 
shaded and hatched bars represent the full sample. One tick mark 
on the y-axis is equivalent to one binary. The classifications CC, 
WW, CW, and WC refer to binaries composed of classical and 
weak-line T Tauri stars with the first and second character de- 
scribing the primary and secondary component. No systems with 
two accreting components were detected. 



4.2.2. Differential disk evolution 

Fig. [T7] shows the separation of Cha I binaries classified by their 
accretion status as CC (both components accreting), WW (no ac- 
creting components), CW and WC (mixed pairs with the primary 
and secondary component accreting, respectively). While the 
majority of binaries in the restricted sample were found to have 
no accreting components (12^^/17 = 70^|^% WW), there are no 
CC binaries (0^^) detected and all accreting components are part 
of mixed systems with either an accreting primary (4^2><CW) or 
secondary (l^^xWC). A qualitatively identical picture is posed 
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by the full sample. The population of these classes agrees well 
with random pairing of the number of individual accretors and 
non-accretors in the restricted sample. About one CC binary 
would be expected, -four mixed pairs, and the rest (-12) are 
expected to be WW. This indicates that the existence of a disk 
around one Cha I binary component has no impact on the prob- 
ability of finding another disk around the other component, and 
both components can be treated as independent in terms of their 
disk probability. The eff^e ct of disk synchroniza tion that was de- 
tected for ONC binaries (iDaemgen et al.l 120121) is not observed 
here. This is probably because Cha I is older than the ONC. 
However, a preference for mixed systems with an accreting pri- 
mary can be observed. WC systems appear more scarce than CW 
as obser ved in previous stud ies of Taurus and other star- forming 
regions (JMonin et al .112007 *). but it was not confirmed by, e.g., 
data of the ONC (4xCW and 3xWC: lDaemgen et al.l[2012h . 

The paucity of WC systems in Cha I points to either a bias 
against their detection or a diff'erent formation and/or evolution 
of secondary disks as a function of the cluster environment. 
Since binaries in this sample were selected from an adaptive 
optics search for companions, no preference for or against ac- 
creting secondary components is expected. The accretion mea- 
surements might even be in favor of the detection of accretion 
in the fainter secondary spectra since the detection limit of the 
mass accretion rate is lower the fainter the system. In fact, the 
only accreting secondary component in this study has the low- 
est detected mass accretion r ate, be low the upper limits of most 
other components (see Sect. 14.2.31 and Fig. [19]). The fact that, 
independently of spectral type, secondaries systematically show 
no or only very weak accretion indicates that we see a real eff'ect 
of the binary environment rather than a selection eff'ect. A real 
paucity of mixed systems with accreting secondaries in Cha I is 
therefore the most plausible explanation. Assuming that almost 
all binaries start in a CC configuration (more than 60% of all 
binaries in the very young 1-2 My r Taurus association are CC 
White & GhelllOOl.: .McCabe et al]l2006l: iKenvon & HartmannI 



19951) . this observation points to a more rapid evolution of the 
disk around the less massive component. 

While the above considerations show that the spectral type 
and thus the absolute mass of the binary components is unlikely 
to cause a preference for mixed systems, we explored the pos- 
sibility that the relative disk lifetimes of primaries and secon- 
daries are governed by the relative masses of the stellar compo- 
nents, i.e., the component mass ratios q = M^Qc/^pnm- Fig.fTSl 
shows the mass ratio as a function of binary separation, indicat- 
ing the accretion type as well. The diagram features mixed sys- 
tems only with mass ratios below ^ = 0.8, while 5 out of 12 WW 
systems measure ^ > 0.8. This suggests a mass ratio-dependent 
mechanism: owing to its smaller Roche lobe size, disk trunca- 
tion is stronger for the less massive component, which may lead 
to a shorter disk li fetime due to an early- ceasing angular mo- 
mentum transport (lArtymowicz & Lubowl 1 1994). A reduction 
of secondary disk lifetime agrees with theoretical expectations 
(it scales with disk radius R roughly as R-^-^-^ for a viscous 
disk), although a significant diff'erence between primary and sec- 
ondary disk lifetime is only expected for mass ratios q < 0.5 
dMonin et al. 2007). 

While mass ratios together with viscous disk evolution can 
explain at least part of the observed CW systems, the more mas- 
sive component is always predicted to keep its disk for longer. 
Accordingly no WC systems would form. That these systems are 
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^^ T 39 B is a spectroscopic binary candidate, which might explain its 
higher inferred mass. 



Fig. 18. Mass ratio q = M^Qc/^prim as a function of binary sep- 
aration. Red open circles show binaries with two non-accreting 
components (WW), blue triangles (up) show mixed systems with 
an accreting primary and green triangles (down) with an accret- 
ing secondary from the restricted sample. Gray symbols repre- 
sent targets that are in the full sample but not in the restricted 
sample. There are two binaries (Hn4, T39 ABPI both WW) 
where the mass of the fainter component is calculated to be 
higher than the primary's mass. The inverse q~^ is plotted for 
these. 



nevertheless observed might be partly due to projection eff'ects. 
Disks in very wide systems that only appear close due to projec- 
tion may evolve without notice of the other star if their outer tidal 
truncation radius is larger than the size of a typical disk, i.e., for 
binary separations »1000au. On one hand, unrelated systems 
can form in any configuration since the lifetime of an individual 
disk can be anywhere between < 1 and ~ 10 Myr. On the other 
hand, WC systems might even be preferred for very wide and 
unrelated components, since it was found that the frequency of 
disks around single stars declines with incr easing stellar mass 
(e.g. jLada eTaDl2006l: [Carpenter et al.ll2006h . 

4.2.3. Mass accretion in Cha I binaries 

Mass accretion rates are shown as a function of stellar mass in 
Fig. [T9j The new Cha I bina ry data are compared \y ith single 
star mass accretion rates from lAntoniucci et aP (1201 ll) . Because 
these were obtained from the same indicator, Bry emission, 
no systematic biases are expected. The distributions of the two 
populations, Cha I binary components and single stars, qualita- 
tively agree in the log(M)-log(M:,) diagram. There are, however, 
too few accretors in the binary sample to statistically confirm 
or reject agreenient wi th the relation M oc M^^^ proposed by 
lAntoniucci et al.1 (l201 ih . 

Interestingly, all accreting components in this sample were 
observed at masses >0.32Mo, i.e., 100% of all accretors are 
among the most massive -60% of all components in the binary 
sample. A KS-test returns a 87% chance that the masses of ac- 
cretors and nonaccretors are drawn from diff'erent parent popu- 
lations. This is neither due to limited sensitivity - the detection 
limits are mostly below the log(M)-log(M*) relation of accretors 
by a bout 0.5 magnitudes - no r a general property of Cha I accre- 
tors: Antoniucci et al.l (1201 ll) did detect single stars with masses 
< 0.32 M© with significant accretion. This corroborates the pre- 
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Fig. 19. Mass accretion rate as a function of stellar mass. Filled 
circles and diamonds are significantly accreting primaries and 
secondaries of Chal binaries. Open symbols are upper limits of 
undetected binary components. Colored symbols denote the re- 
stricted sample, gray symbols the full sample only. Green filled 
triangles show C hal single stars also measured from Bry by 
I Antoniucci et al.ii (201 1). These authors also derived a power-law 
fit to their data, included as a green dashed line. Upper-limit val- 
ues of the I Antoniucci et all non-detections are not published. 

viously inferred low relevance of spectral type or mass for the 
existence of accretion signatures of secondary stars: since low- 
mass accretors < 0.32 M© do exist in Chal singles, but none are 
observed in binaries (neither primary nor secondary), it is not the 
lower mass of the secondary that reduces the chance of dealing 
with an accreting star, but the presence of a more massive stellar 
companion is the most likely cause for the disappearance of the 
disk around the secondary. 

While these considerations suggest a dependence of accre- 
tion parameters on binary parameters, the small number of bina- 
ries in Cha I and the low fraction of accretors limits the statistical 
significance of these results in ChamaeleonI. Future exploration 
of mass accretion rates must therefore rely on the compilation of 
data from various star- forming regions. 

5. Conclusions 

We have presented an observational near-infrared study of 
the individual components of 19 TTauri binary stars and 7 
triples, including one newly resolved tertiary component, in the 
ChamaeleonI star forming region. The target sample covers 
close to all currently known multiples in the separation range 
between 0'.'2 and 6'\ Imaging data in JHK^L bands were taken 
with VLT/NACO together with close-in-time 1 .5-2.5 jam integral 
field spectroscopy with VLT/SINFONI, covering the accretion- 
induced Brackett-y feature. The combined data set enabled us to 
measure the relative position, magnitude, color, spectral type, 
and the equivalent width of Bry emission of each individual 
component. Stellar luminosities and eff'ective temperatures were 
calculated, which enabled us to derive stellar component masses, 
ages, and radii through a comparison with pre-main sequence 
evolutionary tracks. This information was combined with the de- 
rived probabilities for the existence of hot circumstellar dust, as 
inferred from K^-L color excess, and accretion disks from Bry 
emission. The data enable us to conclude the following: 



1. The components of a binary are more coeval than random 
pairs of Chal binary components, with a significance of 
2.3a. 

2. Components of binary stars in Chal show signs of accre- 
tion less often than single stars in the same region. From 
the new data of TTauri binaries with separations of 25- 
1000 au, a fraction of F = 19^^% of all stellar compo- 
nents has been inferred to harbor an accretion disk. This 
is less than half of the single star disk fraction in Cha I of 
^single _ 44 + 8%^ Apparently, the presence of a stellar 
binary component causes accretion to cease earlier. A re- 
duction of the accretion disk frequency is particularly pro- 
nounced in close binary systems of separations 25-100 au 
where the accretion disk frequency i^<iooau = 10^^^% is less 
than half of the accretion fraction in wider binaries (100- 
1000au,F>iooau = 23!ii%). 

3. The dust disk fraction of binary components is indistinguish- 
able from that of single stars, g = 50^^^ % of all binaries 
with separations between 25 and lOOOau were found with 
K^-L excess, identical to the measured value in singles. The 
frequency of hot circumstellar dust around the components 
of binary systems with separations < 25 au, in contrast, is re- 
duced to < 25%. 

4. All accreting components were found to be members of 
mixed pairs of accreting and non-accreting components. The 
presence of a disk around either binary component appears 
to be unrelated with the disk around the other component, 
suggesting no synchronized disk evolution. Mixed systems 
were not found with equal component masses, but their 
mass ratios are always below ^ = 0.8. Pairs without accret- 
ing components occupy a range of mass ratios including 
close-to equal-mass binaries with ^>0.8. Unequal compo- 
nent masses and the induced stronger disk truncation of the 
less massive component are a possible explanation for the 
diff'erential disk evolution. 

5. Out of the seven identified accretors only one is a secondary 
star, i.e., the less massive star in the respective binary sys- 
tem. This is also the weakest accretor found in the sample. 
This suggests that accreting secondaries are rare and/or only 
weakly accreting. The existence of single stars with masses 
comparable to those of the binary components but signifi- 
cantly stronger accretion indicate that this is a consequence 
of the presence of a higher-mass stellar companion to a disk- 
hosting star. 

6. The measured mass accretion rates of binary components 
agree well with the mass accretion rates of single stars 
of Chal for components >O.3M0. Accreting binary com- 
ponents with masses below 0.3 M©, however, are less fre- 
quently observed than around single stars. 
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